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ABSTRACT 

The anisotropy of 1. 3 - 2. 3 MeV protons in interplanetary 
space has been measured using the Caltech Electron/Isotope Spec- 
trometer aboard IMP-7 for 317 6-hour periods from 72/273 to 74/2. 
Periods dominated by prompt solar particle events are not included. 
The convective and diffusive anisotropies are determined from the 
observed anisotropy using concurrent solar wind speed measure- 
ments and observed energy spectra. The diffusive flow of particles 
is found to be typically toward the sun, indicating a positive radial 
gradient in the particle density. This anisotropy is inconsistent with 
previously proposed sources of low-energy proton increases seen at 
1 AU which involve continual solar acceleration. 

The typical properties of this new component of low-energy 

cosmic rays have been determined for this period which is near 

solar minimum. The particles have a median intensity of 0. 06 pro- 
2 

tons/cm -sec-sr-MeV and a mean spectral index of -3. 15. The 
amplitude of the diffusive anisotropy is approximately proportional to 
the solar wind speed. The rate at which particles are diffusing to- 
ward the sun is larger than the rate at which the solar wind is con- 
vening the particles away from the sun. The 20 to 1 proton to 
alpha ratio typical of this new component has been reported by Me- 
waldt, et al. (1975b). 

A propagation model with K rr assumed independent of radius 
and energy is used to show that the anisotropy could be due to in- 
creases similar to those found by McDonald, et al. (1975) a t ~ 3 AU. 
The interplanetary Fermi-acceleration model proposed by Fisk 
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(1976) to explain the increases seen near 3 AU is not consistent with 
the ~ 12 per cent diffusive anisotropy found. 

The dependence of the diffusive anisotropy on various parame- 
ters is shown. A strong dependence of the direction of the diffusive 
anisotropy on the concurrently measured magnetic field direction is 
found, indicating a Hj. less than Hy to be typical for this large data 
set. 
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I. INTRODUCTION 

Observations of ~ 1 MeV protons in interplanetary space 

have revealed rich phenomena with the flux of these particles vary. 

ing over many orders of magnitude and the time scales for these 

changes ranging from hours to days. It is to investigate the source 

of these particles and their variations that the present study has beer 
undertaken. 

McCracken and Rao (1970 ) have reviewed the observed char, 
acteristics of these particles. The authors group individual increas- 

events, into two classes prompt events and delayed events. 
Prompt events, sometimes referred to as "classical- solar particle 
flares, have several distinct characteristics. They have short rise 
times, typically hours, and longer decay times, typically tens of 
hours. Velocity dispersion is observed during the onset of the event: 
the flux of higher velocity particles increases before the flux of low- 
er velocity particles. The onset of the prompt event normally fol- 
lows an enhancement of activity on the sun including optical flares, 
x-ray emission, and radio emission. A large directional dependence 
of the flux is observed early in these events with most of the parti- 
cles streaming along the interplanetary magnetic field from the sun. 
These prompt events are generally thought to be due to nearly im- 
pulsive acceleration and injection of particles into the interplanetary 
medium by solar flares. Detailed studies of the anisotropies during 
prompt events have been done by McCracken, et al. (1971 ), Rao, 
etal. (1971), and Allum, et al. (1974). These studies have indi- 
cated the importance of anisotropy measurements for understanding 
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low-energy proton events. The present work will extend the work to 
periods between prompt events. 

Increases that do not have the characteristics of prompt events 
have been termed delayed events by McCracken and Rao. These 
events typically have a longer rise time than prompt events, and no 
velocity dispersion is observed. Delayed events are also smaller 
than the larger prompt events. 

Some delayed events occur during the decay phase of prompt 
events {Rao, et al, 1967a ; Lin, et al. , 1968 ). The present work is 
restricted to periods between prompt solar events, so such aelayed 
events are not included. 

Other delayed events are not so clearly associated with prompt 
events. Bryant, et al. (1965 ) found events lasting several days that 
recurred every 27 days (one solar rotation period) for 7 successive 
solar rotations. The low-energy threshold was 3 MeV. These 
events were pictured as approximately steady-state streams of par- 
ticles co -rotating with the sun, whose observed time development 
was caused by their rotation past the observer. 

Periods of ~ 1 MeV proton enhancements lasting - 10 days 
have been reported by Fan, et al. (1968 ). Anisotropies were re- 
ported for two periods of ^ 6 hours each, and the streaming was 
found to be coming from the sun's direction. These enhancements 
were interpreted as evidence of continual acceleration, near -sun 
storage in regions extending ~ 180 in solar longitude, and gradual 
injection into the interplanetary medium. 


Anderson (1969 ) reported events not associated with prompt 
events similar to those found by Rao, et al. (1967a) during the decay 
of prompt events. Anderson proposed that these events were co~ 
rotating with the sun and the particles were escaping from a storage 
region near the sun. 

Further analysis of the events reported by Bryant, et al. 
(1965) was done by McDonald and Desai ( 1971 ) , who interpreted the 
events in terms of storage near the sun. An upper limit of 20 per 
cent was put on possible anisotropies. 

Complex variations in the flux of protons at energies 2:0. 3 
MeV were reported by Krimigis, et al. (1971 ) with large (/t 50 per 
cent) anisotropies, indicating essentially continuous flow of particles 
from the sun. Roelof and Krimigis (1973 ) modelled these observa- 
tions as scatter -free propagation from the sun to the earth with the 
variations caused by either connection to different regions on the sun 
or time variations in the rate of injection. 

Although the delayed events discussed above have varied time 
developments, the models used to explain them are similar in many 
ways. Instead of the nearly impulsive acceleration and injection 
found in prompt events, the delayed events have been pictured as re- 
sulting from either continuous acceleration and injection or intermit- 
tent acceleration and storage between the sun and the earth. Thus, 
the sun is viewed as a continuing source of low-energy protons even 
during periods when prompt events are not observed. Anisotropies 
consistent with this solar source have been measured for some 
events, as indicated above. 
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A more comprehensive study that also suggested continuous 
acceleration was reported by Kinsey (1970 ) , who examined proton 
fluxes in the energy interval 4 - 80 MeV during the time interval 24 
May 1967 to 20 August 1968 and found a highly variable component at 
low energies. He suggested this component is due to a continuous, 
but variable, solar source. Thus, the sun is seen as a continuous 
source of particles whose flux is occasionally large enough to be ob- 
served as individual events. 

The above observations present a rather consistent view, but 
there are two observations which do not fit. Rao, et al. (1967b), re. 
ported a nearly zero anisotropy for protons in the energy interval 
7. 5 - 45 MeV during extended periods of 1965 and 1966 when no 
events were observed. As discussed by Forman and Gleeson (1975) , 

this result is difficult to interpret using contemporary propagation 
theory. 

A recent observation by McDonald, et al. (1975 ) indicates the 
existence of at least some periods when flow back toward the sun 
might be expected. They report that larger delayed events were seen 
by Pioneer 11 between 2 and 5 AD than seen at earth during the same 
period. These events were interpreted as co-rotating streams popu- 
lated by interplanetary acceleration of energetic particles. 

The present study will report the observed anisotropy of 1. 3 to 
2. 3 MeV protons during 317 6-hour periods for a time span from 
72/273 to 74/2, omitting periods when prompt events are observed. 

The data will be interpreted in light of simultaneous measurements of 
the solar wind speed and the interplanetary magnetic field direction 
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using the propagation theory presented in Chapter II. These observa- 
tions directly test the models discussed above that propose that the 
enhancements of the low-energy proton flux seen at 1 AU are caused 
by quasi -continuous injection of particles by the sun. 

The present study differs from those discussed above in com- 
bining the following features. First, many periods from a long time 
span are used. Previous studies have typically been of a few events 
selected by the investigator. Second, for each period the anisotropy 
due to the effect of the solar wind has been subtracted from the ob- 
served anisotropy, leaving the anisotropy due to diffusion. For the 
first time, the dependence of this diffusive anisotropy on such pa- 
rameters as the magnetic field direction is directly determined. 

Third, the background of the Caltech instrument makes a negligible 
contribution to the results reported in this work. Figure 1-1 com- 
pares the fluxes used in the works cited above to those used in the 
present study. Because the observations have been taken at some- 
what different energies, the differential flux at 1 MeV indicated by 
the observations hat been used. An energy spectrum of dj/dT = kT’ 3 , 
where j is the differential intensity and T is the kinetic energy, 
typical of spectra observed at energies near 1 MeV, has been assumed 
for calculating the differential flux. Some of the previous studies 
were at fluxes much higher than those used in the current study. In 
addition, the quiet-time anisotropy measurement by Rao, et al. 

(1967b) is seen to have a higher equivalent flux than some of the pe- 
riods used in this work. Thus, the current measurements are rele- 
vant to understanding this previous result. 
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Figure 1-1 

A comparison of proton fluxes for delayed events and quiet times for 
various experiments. The references for the designated experiments 
are: 

UTD-1 Rao, et al., (1967a). 

UCB-1 Lin, et al., (1968). 

GSFC-1 Bryant, et al., (1965); McDonald and Desai (1971). 
Chicago Fan, et al., (1968). 

UCB-2 Anderson, (1969). 

JHU/APL Krimigis, et al., (1971). 

GSFC-2 Kinsey, (1970). 

UTD-2 Rao, et al . (1967b). 

Caltech present work. 
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II. EXPERIMENT 

A. Spac e craft 

The energetic particle data used in this study were obtained 
from the Caltech Electron/Isotope Spectrometer (EIS) experiment 
aboard the IMP -7 spacecraft. IMP -7 was launched into an orbit rang- 
ing from 32 to 36 earth radii in September, 1972. The orbit is in- 
clined to the ecliptic by 28°. A rotation of this orbit into the ecliptic 
plane is shown in Figure II- 1, along with the average position of the 
earth's magnetosphere (Behannon, 1968 ). The satellite is sunward 
of the earth’; 3 bow shock during the majority of its orbit. The space- 
craft is spin stabilized with a rotation rate of ~45 rpm. The spin axis 
is within 2° of the South Ecliptic Pole; the Caltech EIS experiment is 
mounted so that it scans the ecliptic plane as the satellite rotates. 
Signals are generated by the spacecraft to indicate the current orien- 
tation of the satellite. Each rotation is divided into eight equal sectors 
as shown in Figure II-2. The signal for Sector 0 is initiated by the 
spacecraft's sun sensor detecting the sun. Because of the position of 
the Caltech experiment relative to the sun sensor, the Caltech experi- 
ment is not viewing the sun during Sector 0 . The signals for Sectors 
1 through 7 are determined by counting clock pulses after the sun- 
sensor pulse. The rate of the clock pulses is adjusted by the space- 
craft so that the 8 sector signals are of equal length and include the 
entire rotation. Tests of the accuracy of the sectoring system are 
described in Chapter IV. 

B. Instrument 


A brief description of the aspects of the Caltech experiment 


{ 
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Figure II -1 

The orbit of lMP-7 rotated into the ecliptic plane. Also shown are 
the average positions of the earth's bow shock and magnetopause as 
determined by Behannon (1968). 



Figure II-2 


Diagram of the IMP-7 satellite as seen from the North Ecliptic Pole. 
The relative positions of the spacecraft sun sensor and the Caltech 

Electron/ Isotope Spectrometer are shown. The average viewing angle is 
given for each of the 8 sectors. 
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IMP-7 SECTORING 


Sector 4 



Solar Ecliptic Coordinate System 
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relevant to the current study is given here; more complete descrip- 
tions are available elsewhere (Hurford, et al. , 1974; Mewaldt, et al. , 
1975a). A cross section of the Caltech EIS telescope on IMP-7 is 
shown in Figure II-3. The telescope consists of 11 fully-depleted 
silicon surface-barrier solid-state detectors and an anti -coincidence 
scintillator viewed by a photomultiplier tube. Each of the solid-state 
detectors except DIO has a pulse height analyzer and discriminator; 
DIO has only a discriminator. The solid-state detectors are all 
nominally 1000 microns thick with the exception of the 47 micron 
thick D2. The telescope is covered with a 2.4 mg/cm* - aluminized 
mylar window. 

The instrument has several modes of analysis. The data used 

in this work come from the narrow geometry mode in which detectors 

DO, Dl, D3, D4, DIO, and Dll are in anti -coincidence. An example 

of a narrow geometry signature is an event for which only D2 triggers. 

This work is concerned primarily with events of this signature, 

2 

which have a geometric factor of 0.21 ±0.01 cm -sr and an opening 
half-angle of 29° ± 1° (Hartman, 1973 ). 

Due to the limited telemetry rate available to the experiment, 
only part of the information produced by the particles which trigger 
the telescope is transmitted to earth and thus recorded. A more de- 
tailed description of the Caltech experiment is available elsewhere 
(Garrard, 1974 ). There are two types of transmitted information of 
interest to this work -- rate information and analyzed event informa- 
tion. 

The rate information consists of the rates at which individual 


! 
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Cross Section of the 


Figure II -3 

IMP-7 Electron/Isotope Spectrometer Telescope. 


> 
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detectors trigger and selected combinations of detectors trigger. 
Table II- 1 lists the available combinations and their mnemonics. 

The rates are deteimined by counting events for a known time inter- 
val; this accumulation is done in the satellite external to the Caltech 
experiment. The accumulation system and the Caltech experiment 
are connected by 9 rate lines labelled A, B, C, E, FI, F2, F3, F4, 
and F5. There are four states of this system, called s ub commutation 
states, each lasting 20.48 seconds. During each state a selected set 
of 9 rates is accumulated. Table II-2 indicates which rates are ac- 
cumulated during the four subcommutation states. The A, B, C, and 
E rates are accumulated separately for each of the 8 sectors using 4 
sets of 8 accumulators. These rates are accumulated for only 14 
spacecraft rotations (~ 18 seconds) out of the 20.48 seconds of the 
subcommutation state to insure that only complete rotations are used. 
The rotation period is measured so that the accumulation time is a 
known, but slowly varying, quantity. Because events are accumulated 
by sectors, anisotropies can be computed for these rates. An exam- 
ple of a sectored rate is PLO, which uses rate line A and is accumu- 
lated during subcommutation states 1 and 3. The PX^D rate shares 
the accumulators associated with rate line A with the ELO rate. The 
remaining rate lines FI, F2, F3, F4, and F5 are not accumulated by 
sectors, and so require a totcil of 5 accumulators. Rates using these 
lines are accumulated continuously. The rates associated with lines 
FI and F2 are read out twice every subcommutation state; the rates 
associated with lines F3, F4, and F5 are read out four times every 


subcommutation state. 


Table II-l. Rates 


Rate 

Logic Requirement 

Nominal Physical Significance 

ELO 

N # S*d2TT*D5-d5h*35'»3T*aip 

Electrons, ^0.2 to *1 MeV, including Compton recoils 

EHI 

N D5*3?F*D6 

Electrons, ^1 to ^3 MeV, including Compton recoils 

PLO 

N*S*D2H*HF«3^*3T‘aTp 

Nuclei, 1.2 to 2.4 MeV + 

PHI 

N*S*D2.D5H*cRr 

Nuclei, 4 to 13 MeV 

DO* 

DO-S 

Electrons, ^0.16 to ^5 MeV, plus nuclei, 1 to 43 MeV 

DOT* 

D0*D1 «S 

Electrons, ^1 to ^5 MeV, plus nuclei, 13 to 43 MeV 

PEN 

N-D5-D6-D10-HTT 

Electrons and nuclei that penetrate the telescope 
(electrons > 3 MeV, nuclei > 30 MeV) 

NEUT 

N*S-d!T-D7 

Neutral particles, such as y-rays, whose Compton recoil 
electrons are detected in D7 

ADC 

HAZ 

D0,D1,,,D11 

D0+D1+...+D9 

The logical OR of all pulse-height-analyzed detectors 
The rate at which the HAZARD flag is set 
Singles rates for individual detectors 


^Nuclei energy limits are for protons. 


n = <iD«<rr»cnr*<Tr 

s - an^arr 

AIP = analysis in progress 


Table II-2 . IMP-7 Rate Block 


A 

B 

C 

ELO 

PHI 

DO 

PLO 

EHI 

DO* 

ELO 

PHI 

D5 

PLO 

EHI 

D2H 


E 

F]_ 

F2 

ADC 

D1 

D3 

DOT* 

D6 

07 

DU 

D8 

D9 

PEN 

HAZ 

D2 


£3 £4 £5 

04 NEUT DO 
DIO NEUT 00* 
D5H NEUT 05 
010 NEUT D2H 
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The other type of information transmitted consists of analyzed 
events. The information for each analyzed event includes which de- 
tectors triggered, the sector during which the event occurred, and 
two pulse heights giving energy -loss information. The pulse heights 
transmitted depend on the signature of the events; the D2 and D5 
pulse heights are transmitted for D25 events, for example. The data 
from one analyzed event are transmitted every 0. 64 seconds. The 
selection of events transmitted is determined by a 5 -level priority 
system. The equations determining an event's priority are given in 
Table II-3. Only events of the two highest priorities, PO and PI, are 
used in this study. These events nominally consist of narrow geome- 
try electrons and nuclei. If the previous event read out was a narrow 
geometry electron, the highest priority is assigned to narrow geome- 
try nucleons. Otherwise, the highest priority is assigned to narrow 
geometry electrons. 

Thus, both the D2 analyzed events and the PLO rate are a 
sample of D2 events. The live time for the PLO rate is ~ 50 percent 
of the elapsed time because the PLO rate is accumulated during 2 of 
the 4 sub commutation states. The live time for D2 analyzed events 
depends on the rates of the different types of narrow geometry events. 
Typically at PLO rates < 0. 5 /second, there are more D2 analyzed 
events than events counted for the PLO rate, while the converse it? 
true at higher rates. 
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P0 = 
PI = 
P2 = 
P3 = 
P4 = 
where 


Jable II-3. Priority Levels 


N«S*c!F* {35F»D5*3?FTe + (D2H + D5H 
N - S * 3TT* {d2h*D5*cT5h«RE + (D2H + D5H 
S*dh* (DO + N*3T'D7) 

N*3F* (D2H + D5) -Dll 
N* (D2H + D5) -CH 


)«RE> + N*3TT.3IT.D5*D10*RE 
)' re> + W*dTT*Hh‘D5*D10‘re 


N = H0.aT-3T.3? 

s = m-m 

OH = hazard bit 

RE = recent electron bit 
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III. COSMIC-RAY PROPAGATION 
A. Diffusion - Convection Model 

The basic principles of cosmic -ray propagation have been re- 
viewed by Jokipii (1971). Propagation is controlled by the interplan- 
etary medium which is a highly conductive plasma expanding nearly 
radially outward from the sun. This solar wind has a bulk velocity of 
~ 300 - 600 km/sec and a flow direction within ~ 5° of radial. Im- 
bedded in the plasma is the interplanetary magnetic field which is 
carried out from the sun by the solar wind. Due to the sun's rota- 
tion, the average field direction forms an Archimedian spiral. The 

average field direction at 1 AU is about 45° from radial. The magni- 

-5 i 

tude of the magnetic field is ~ 5 XlO gauss at 1 AU; a l^r MeV proton 

9 -4 

has a Larmor radius of 3.4 X 10 cm or 2. 3 X 10 AU in this field. 

Individual measurements of the interplanetary magnetic field 
fluctuate from the mean value. These fluctuations cause the cosmic 
rays to scatter from simple helical motion. After many such scat- 
terings, propagation can be described as diffusion with a diffusion 
tensor x • Because of the presence of the magnetic field, diffusion is 
in general not isotropic. Assuming that there is no difference in the 
two directions perpendicular to the field, the diffusion tensor has the 
form 

( • *i - h t 

ht hi 
0 0 

where the z-axis is taken to be parallel to the direction of the mag- 
netic field. K|| and Hj. characterize diffusion parallel to and per- 
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pendicular to the magnetic field, respectively. Hr describes the 
streaming produced by density gradients perpendicular to the mag- 
netic field. 

No generally accepted theory exists for computing H from 
observed properties of the interplanetary magnetic field. The quasi- 
linear approach used by Jokipii (1966) to determine x from observ- 
sble statistical properties of the magnetic field requires certain sim- 
plifying assumptions whose validity has been questioned (Fisk, et al. , 
1974; Birmingham and Jones, 1975). The results do suggest that the 
important quantity is the amount of power in the fluctuations of the 
magnetic field at wavelengths comparable to the Larmor radius of the 
particle being scattered. This resonant scattering theory indicates 
that Hx would be much smaller than H|j if it weren't for an additional 
term contributing to x^ that has been interpreted as being due to the 
random walk of the field lines (Jokipii and Parker, 1969). That is, 
neighboring field lines at the sun can be separated at 1 AU allowing, 
for example, solar flare particles to have a larger azimuthal extent 
at 1 AU than at the sun without any scattering perpendicular to the 
field lines taking place. The relative size of x x and K| ( remains un- 
resolved. Values of K x / K\\ ranging from 0 (Wibberenz, 1974) to - 1 
(Jokipii, 1971) have been suggested. 

B. Anisotropies 

The following picture of particle propagation emerges. The - 1 
MeV cosmic-ray protons execute helical motion about the average 
magnetic field as they are being convected outward from the sun with 
the field at the bulk velocity of the solar wind. In addition, the cosmic 
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rays are scattered by the irregularities in the magnetic field. For. 

man and Gleeson (1975, have shown that the differential streaming is 
given by 

S = CVU - VU , ^ 

where U is the differential number density, V is the solar wind ve- 
locity, and C is the Compton - Getting factor 


C = <2-Oy)/3 , 


(3-2) 


where a = (T + 2m c 2 ,/(T + m c 2 , and y is the spectral index given by 

v = j/a&iT 3) 

where j is the differential intensity, j = wU/4tt , and w is the par- 
ticle velocity. For the low energy protons relevant to this study, 
a —2 . The vector anisotropy is dimensionless and indicates the reU 
tive amount of streaming toward different directions. It is defined b 
= 3S-/UW = | [CV-^vUJ/U], r cON+ r DIF . (3-4) 

The anisotropy is the sum of two terms, a convective term 
and a diffusive term. The convective term is proportional to the 
solar wind velocity and is related to the energy spectrum of particles 
through the Compton - Getting factor. The diffusive term is deter- 
mined by the product of the diffusion tensor and spatial gradients in 
the particle density. Thus, to learn about the spatial distribution, 
the convective term must be subtracted from the observed anisotropy. 
The typical size of the convection term is ~ 20 per cent. 

The diffusion - convection model has been used by McCracken, 
et al. (1971) to explain the observed time development of the anisotro- 
py during prompt solar particle events. As discussed in Chapter I, 
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Figure III- 1 ! 

The model for the evolution of the anisotropy during a prompt solar 
particle event. Early times are periods 5 1 day after the onset of the 
event, late times are periods from ~1 to ~4 days after then onset, and 
very late times are periods >~4 days after onset. The figure is adapted 

from McCracken, et al., (1971). 
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study. All the quantities affecting r cON are known. This means 
£ dif can be determined by subtracting from . The 

quantities in £ D j f are not as well known. Wibberenz (1974) has 
summarized the estimates of * rr determined from observed time, 
to-maximum-intensity for prompt solar events. A typical value for 

~ MeV P rotons is 3 x 1° 2 ° cm 2 /sec. Evidence of gradients some- 
times as large as + 300 per cent/AU has been reported by McDonald, 
et al. (1975). Substituting these values into eq. (3-4) indicates dif- 
fusive anisotropies of ~ 10 per cent toward the sun could be observed 
at least during certain periods. 
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IV. DATA ANALYSIS 

A. Particle Selection 

The anisotropy measurements are from the Caltech EIS ex- 
periment described in Chapter II. The anisotropies are computed 
using analyzed events which trigger only D2, which will be referred 
to as D2 analyzed events. Anisotropies using the FLO rate are used 
as a consistency check. The energy intervals of incident protons which 
produce the observations are: 

_ Nominal 73/86 - 73/154 

D2 analyzed events 1. 33 to 2 . 32 MeV 1. 24 to 2. 17 MeV 

PLO 1. 18 to 2.37 MeV 1. 12 to 2. 37 MeV 

A change of ~ 4 channels in the offset of the D2 pulse height analyzer 
lowered the n omin al energies during the period 73/86 to 73 / 

154. The upper limit for energy loss for D2 analyzed events is cho- 
sen to eliminate particles that lose too much energy in D2 to be pro- 
tons. This limit eliminates ~ 5 per cent of D2 events. An upper limit 
to the electron contribution to D2 events can be determined using the 
method presented by Lupton and Stone (1972) for determining electron 
detection efficiencies in solid-state detectors. The authors plot the 
maximum efficiency, that is, the efficiency at the incident electron 
energy at which the efficiency is highest, as a function of the energy 
threshold of a 50 -micron detector. Using the energy threshold for D2 
analyzed events, the maximum efficiency is < 10“ 4 at all times. The 
flux of electrons is monitored by the ELO rate which has a geometric 
factor about one ^fourth that for D2 events at incident electron energies 
~ ^ MeV. The ELO rate is ~ 0, 15/second for the periods used in this 
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c -5 

study, so the electron contribution to D2 analyzed events is ^ 4X 10 / 

second, or < 1 event in 6 hours. Thus, electron contributions to D 2 
analyzed events have essentially no effect on the computed anisotro- 
pies. 

B. Anisotropy Determination 

Anisotropies are computed from 6-hour averages of D2 ana- 
lyzed events during the time interval 72/273 to 74/2. Six-hour aver- 
ages are used to accumulate enough events so that statistical uncer- 
tainties will not be large compared to the measured anisotropy for a 
large fraction of the available periods. As shown in Figure IV- 1, a 
typical PLO rate is ^ 0. 02 /s econd, which produces about 400 D2 
analyzed events in 6 hours. An anisotropy measurement using 400 
counts will have a statistical uncertainty wf ~ 7 per cent compared to 
a typical observed anisotropy of ~ 10 per cent. The method used for 
calculating anisotropies and their statistical uncertainties is presented 
in detail in Appendix A; a brief account is given here. 

The anisotropy is calculated by fitting in a least squares sense 
the observed number of counts by sector to the function: 

f(4>) = A + B cos + C sin (4-1) 

The parameters of the fit A , B , and C are readily solved from line- 
ar equations. Anisotropies usually are expressed by an equivalent 
function, but one with different parameters: 

f(4>) = A q (1 + i cos(<},-<}> 0 )) . (4-2) 

A direct fit to this function leads to non-linear equations for the pa- 
rameters. To avoid this, the parameters of eq. (4-2) are determined 
from the parameters fit to eq. (4-1). £ is then multiplied by 1. 079 to 


Figure IV- 1 


The distribution of the 6 -hour averages of the FLO rate observed by 
the Caltech EIS experiment aboard IMP-7. The flux of 1. 2 to 2. 4 MeV 
protons is computed by approximating the FLO rate as being due only 
to protons stopping in D2. Heavier nuclei and protons penetrating D2 
typically produce - 15# of the PLO rate. 
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correct the smoothing effect due to the finite opening angle of the tele- 
scope and the finite number of sectors. Error bars on £ and 6 
used in this work are defined by 


§ 


I 


1.079 a 
~~2 7 T • 

o 


where flf is average uncertainty in counts in a sector, 


a 





57.296° . 


(4.3) 


(4-4) 


The determination of the probability distribution for the true 
value of £ and <J> o given an observation is discussed in Appendix A. 
Let 


z = %/ts . 

For z £ 2, that is, for statistically significant observations, the 

probability distribution for the true £ is approximately Gaussian, 

centered near the observed value with a standard deviation given by 

eq. (4-3). For z £ , the most likely value of £ is 0 , and eq. 

(4-3) is only an estimate of the statistical uncertainty. 

The observed is always the most likely value for the true 

4> 0 • The C , given in eq. {4-4) is the standard deviation of the prob- 

ability distribution for the true in the limit of large z. For z ^ 

2, ct is only an indication of the uncertainty in the determination of 
v o 

V 


For | « 1 , eq. (4-2) describes a nearly isotropic distribu- 
tion with a small anisotropy of amplitude £ and a maximum flux ob- 
served when the telescope is pointing in the direction <b . The 

o 


streaming direction 



is defined by 
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<b f. - 4> + 180 . 

r £, 

The functional dependence given by eq. (4-2) is found to pro- 
vide a good fit to nearly all the data used in this study. Define a 

2 

goodness -of-fit parameter x, using the data points y. , their sta- 
tistical uncertainties , and the fitting function y(x^) in the usual 
way: 

X* = (y. - y(x.)) 2 /a. 2 ] /v , 

2 

where v is the number oi degrees of freedom. X v is computed for 

each period. If the deviations of the observations from the fit are 

2 

due only to statistical fluctuations, the mean X^ equals 1 . Figure 

2 

IV -2 presents the observed distribution of X v and the distribution 
predicted due to statistical fluctuations. Periods when the PLO rate 
is larger than 0.3/second or smaller than 0.01/second, periods 
when the satellite is not on the sunward side of the earth, and periods 
near the beginning of prompt solar particle events are not included in 
the distribution. These periods are also not included in the final 
data set used in this study; the selection criteria for the final data set 
are discussed in Chapter V. 

Both the observed and predicted distributions peak near 0.7 

and fall to half maximum near 1.4, indicating eq. (4-2) provides a 

good fit to most of the periods. There are more periods with large 
2 

X than predicted, but they comprise a small fraction of the total 

2 

number of periods; only 3 per cent of the periods have x^ a 4.0. 

2 

The mean x of 1. 28 indicates that the average deviation from the 
fit is 1.13 times as large as that due to statistical fluctuations. 

Thus, deviations from the function fit to the observations are domi- 
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Figure IV-2 

The distribution of the goodness-ofrfit parameter determined from 
the anisotropy fits to D2 analyzed events. The mean of this distribu- 
is indicated. The number of degrees of freedom v for the fit is 5. 

The data set consists of periods from 72/273 to 74/2 when the average 
FLO rate is between 0.01/second and 0.3/second. In addition, the av- 
erage longitude of the satellite in Geocentric Solar Ecliptic coordi- 
nates, $ SAT , is required to be between -90° and + 90°. Only periods 
for which the magnetic field B, the solar wind speed V, and spectral 
index y are known are included in the data set. As discussed in the 
text, periods near the onset of identified prompt solar particle events 
are not included. Also shown Is the distribution of expected if the 

deviations from the anisotropy fits were due only to statistical fluctu- 
ations. 
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nated by statistical fluctuations due to the finite number of counts. 

Tests have been made to determine that the sectoring system 
aboard IMP -7 is not malfunctioning and thereby introducing significant 
instrumental errors into the measured anisotropy. The first test 
uses background rates which are known to be nearly isotropic. As 
discussed by Hurford, et al. (1974), the ELO and EHI rates have a 
substantial background rate caused by gamma rays Compton scatter- 
ing in the telescope and triggering D5 to simulate ELO events or D5 
and D6 to simulate EHI events. During most of the time, this back- 
ground rate is larger than the true electron rate as determined from 
DO events. The gamma rays are produced by the interaction of rela- 
tivistic cosmic-rav nucleons in the spacecraft. This process is 
largely independent of the spacecraft orientation, so the resultant 
ELO or EHI events are nearly isotropic. Consequently, if the instru- 
ment is working correctly, the anisotropy computed for these back- 
ground rates should be nearly zero. Anisotropies are computed using 
data from 72/273 to 75/18, rejecting periods when the satellite is not 
sunward of "he earth. Figure IV-3 has histograms of the ELO and 
EHI rates for this period. In order to avoid rates above background, 
ELO anisotropies include only periods when the ELO rate is less than 
0. 053 /second, and EHI anisotropies include only periods when the 
EHI rate is less than 0. 025/ second. Periods near the beginning of 
prompt solar particle events are also excluded. The computed ELO 
anisotropy amplitude is consistent with 0 and is less than 0. 67 per 
cent at a 95 per cent confidence level. The EHI anisotropy is also 
consistent with 0 and is less than 0. 86 per cent at the 95 per cent 


37 


Figure IV-3 

The distribution of the observed ELO and EHI rates. The arrows indi- 
cate the respective upper limits of rates for periods included in the 
calculation of the anisotropies of the backgrounds of ELO and EHI. 
There are 1599 periods in each panel. 
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confidence level. Thus, any long-term bias affects these anisotropies 
by ~ 1 per cent. As discussed in Chapter II, the rate accumulation 
system is shared by different rates through the use of 9 rate lines 
and 4 subcommutation states. Since ELO and PLO share the same 
rate line and EHI and PLO are accumulated during the same subcom- 
mutation states, any instrumental error in the PLO anisotropy should 
also be ~ 1 per cent. 

Although no large instrumental errors were found in the 
anisotropies using rate data, the instrument might introduce errors 
into the anisotropies of analyzed events. To check this, the second 
test compares the anisotropy for DZ analyzed events to the anisotropy 
for the PLO rate. As explained in Section A, these two event types 
are samples of D2 events with slightly different energy intervals. 
Thus, the two event types should measure approximately the same 
anisotropy if the instrument is functioning properly (except at high 
rates -- see Chapter V). The difference vector, 0 ” 

Jr -LjVJ , Dii? 

*■# 

^D2 DIF’ * S P^ otte< * Figure IV-4 for 6-hour periods when the PLO 
rate is between 0.03/second and 0. 3 second, and the satellite is sun- 
ward of the earth. Periods near the beginning of prompt solar events 
have not been included. The mean differences in the x- and y-direc- 
tions are -0. 3 per cent and 0. 5 per cent, respectively, with a sta- 
tistical uncertainty of ± 0. 5 per cent. The average particle velocity 
for DZ analyzed events is ~ 1.04 times larger than the average ve- 
locity for PLO events. Thus, depending on how the spatial gradients 
and the diffusion coefficient vary with energy, there may be a sys- 
tematic ratio between the anisotropies of ~ 1.04; for a typical diffu- 


Figure IV-4 

The difference in t DIF determined using the sectored PLO rate and 
^ PI p determined using D2 analyzed events. Each dot indicates one 
6-hour period. A typical 1 a error bar is indicated. 
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sive anisotropy of ~ 15 per cent, this is a systematic error of ~ 0, 6 
per cent -- comparable to the statistical uncertainties. The small 
mean differences found are consistent with no instrumental bias be- 
tween the anisotropies of analyzed events and rates. 

The two tests together indicate instrumental errors in the 
anisotropy of DZ analyzed events are & 1 per cent. 

As discussed above, the ETO rate is dominated by an iso- 
tropic background much of the time, producing a measured anisotropy 
insensitive to the true electron anisotropy. In contrast, DZ events 
are comparatively free of such background. During days 66 to 70, 
1973, the average PLO rate was 3.5 ± 0. 4 X 10~^/second, placing an 
upper limit on the background rate. This rate is ~ 1 count/sector/ 
6-hour period. The typical PLO rate limits for periods used in this 
study are 0.01/second to 0. 3/second. If the lowest PLO rate ob- 
served were due entirely to background, a real 10 per cent anisotropy 
at a PLO rate of 0.01/second would be reduced to 9.7 per cent. Thus, 
any background contribution to D2 anisotropies will have at most a 
minimal effect. 

C. Energy Spectrum 

As discussed in Chapter III, part of the observed anisotropy is 
understood to be due to the Compton-Getting effect, whose size de- 
pends on the spectral index y given by 

y = d 071 j/ d 071 T 

where j is the differential intensity, and T is the particle kinetic 
energy. The spectral index is determined from data from the Caltech 
EIS experiment. The function 
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dj/dT = AT Y (4-5 ) 

is fit in a least squares sense to the deduced omni-directional inci- 
dent particle spectrum. A detailed description of this procedure is 
presented in Appendix B; a brief summary follows. 

The observed energy losses in D2 due to stopping protons are 
grouped into 9 energy bins. The number of D2 events with these 
energy losses are summed over a 6-hour period. Not all of the en- 
ergy losses observed in these energy bins are due to stopping pro- 
tons. Stopping heavier nuclei, mostly alpha particles, make a small 
(~2-| per cent) contribution. A correction for these heavy nuclei is 
made by extrapolating their observed spectrum at higher energies. 

Some protons that do not stop in D2 also contribute to the D2 
spectrum. Detectors 3, 4, and 5 are used to reject such protons, 
but a thin dead layer around the inside of D3 and D4 allows some pro- 
tons which penetrate D2 to stop in but not trigger D3 or D4. The 

number of such protons is given by: 

T 

2 

N = J F(T)[1 -e (T)1 dT , 

T 

1 

where T is the proton's incident energy, protons with incident ener- 
gies between T^ and T^ lose the appropriate amount of energy in D2, 
F is the fluence of protons, and e is the efficiency with which D3, 

D4, and D5 detect penetrating protons. T^ and T^ are determined 
from range-energy tables, F(T) is determined from observed D25 
events, and e(T) is estimated from results of an accelerator run 
and the instrument's response to large solar particle events. Pene- 
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trating protons typically contribute ~ 10 per cent to the total counts 
observed in D2. 

After the corrections due to heavy nuclei and penetrating pro- 
tons are made, the observed energy loss bins are corrected for the 
energy loss in the mylar window covering the telescope to determine 
the corresponding incident energy bins. Finally, the corrected counts 
and the incident energy bins are fit to the power law energy depend- 
ence given by eq. (4-5). The appropriateness of this dependence is 

investigated using the previously defined goodnes s -of-fit parameter 

2 2 
X . Figure IV-5 presents the distribution of observed X^ and that 

predicted for 7 degrees of freedom. Periods when PLO is larger 

than 0. 3/second or smaller than 0.01/second, when the satellite is 

* 

not sunward of the earth, and periods near the beginning of a prompt 
solar particle event are not included. The closeness of the observed 
mean of 1. 28 to 1. indicates that statistical fluctuations dominate 
observed deviations from a power law energy dependence. 

Figure IV-6 shows the distribution of calculated y's. The 
values range from -1 to -5 around a mean of -3. 15. The finite num- 
ber of counts used to calculate y produces a statistical uncertainty in 
the result which produces a statistical uncertainty in the determina- 
tion of the convective anisotropy; = — ( 1 -y ) V , so 


a £, CON V 


(2V/w)a 

V 


where V is the solar wind speed and w the average particle velocity. 
This uncertainty is typically 1/3 of the statistical uncertainty in the 
determination of the observed anisotropy using D2 analyzed events 
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Figure IV-5 

The distribution of the goodness -of -fit parameter determined 
from the energy-spectrum fits. The mean of the distribution is 
indicated. Also shown is the distribution expected if deviations 
from the fits were due only to statistical fluctuations. 
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Figure IV "6 

The distribution of the spectral index y. Also shown are the mean 
and standard deviation of the ensemble as well as a typical 1 a 
statistical uncertainty in the determination of an individual y. 
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and so contributes only ~ 10 per cent of the uncertainty in the determi- 

nation of the diffusive anisotropy. 

The average particle velocity used in anisotropy equations 

such as (3-4) is determined by: 

t 2 t t 2 

( 1 ) = [J T^T-MT]/[ r T^dT] , 

where T and T ? define the energy interval contributing to the aniso- 
tropy. The resultant average particle velocity is rather insensitive 
to the value of y ; a change of y from -3 to -4 changes > by 
1 per cent for D2. analyzed events. 

D. Solar Wind 

As indicated by eq. (3-4), the solar wind speed and direction 
are needed to compute the convective anisotropy. Hourly averages of 
the solar wind speed measurements by the MIT plasma experiment on 
IMP -7 are combined into 6-hour averages for use in this study. The 
distribution of 6-hour averages is shown in Figure IV -7. The standard 
deviation in the measurements for each of these averages is also calcu 
lated. Figure IV -8 shows the distribution of the ratio of this standard 
deviation to the average. For 90 per cent of the periods the standard 
deviation is less than 5 per cent of the average for the 6-hour period, 
indicating that the average is a good approximation to the solar wind 

speed for the entire 6 -hour period. 

Solar wind speeds larger than 700 km/second have been indi- 
cated as possibly unreliable in the preliminary data set available. 
Consequently, the 134 6-hour periods having an hourly average speed 
greater than 700 km/second are omitted from this study. 


50 


Figure IV -7 

The distribution of the 6-hour averages of the observed solar wind 
speed. The mean 6-hour average is indicated. 
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Figure IV-8 

The distribution of the ratio of the standard deviation and the mean 
of the ensemble of individual solar wind speed measurements used in 
computing the 6-hour average solar wind speeds. 
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The solar wind is assumed to be radial. Previous measure- 
ments of the azimuthal flow angle have been reviewed by Wolfe (1972). 
The mean direction found by different experiments varies from +3.0° 
to -2. 52 from radial, although this variation may be due to system- 
atic errors. A typical standard deviation in the observations of a 
single experiment is 3°. For a typical convective anisotropy ampli- 
tude of -20 per cent, a 3° error in the solar wind direction produces 
an error of ~ 1 per cent in the direction perpendicular to radial and a 
negligible error in the radial direction. 

E. Magnetic Field 

As discussed in chapter III, the diffusion of low-energy pro- 
tons in interplanetary space is understood to be controlled by the mag- 
netic field. Hourly averages of the interplanetary magnetic field have 
been obtained from the National Space Science Data Center for the 
period 72/273 to 74/15. Thu primary source of these data is the Im- 


perial College magnetometers aboard the earth-orbiting HEOS-1 and 
HEOS-2 satellites. The remainder of the data are from the GSFC 
magnetometer aboard IMP-8. 

The field direction used for a 6-hour period is determined by 
computing 

B* = < — L~ hour > 

6-hour • — ' 


B 


1 -hour 


The direction defined by the projection of this field onto the ecliptic 
plane is used as the 6-hour average field direction. The distribution 
of 6-hour average field directions is shown in Figure IV-9. 

Significant variations are found within the 6-hour periods. A 
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Figure IV-9 

The distribution of the 6-hour averages of the observed interplanetary 
magnetic field direction. 
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A 4> is defined for each hourly average field direction: 

“ ^l~hour “ ^6-hour 1 

where the 1-hour period is included within the 6-hour period. Figure 

IV- 10 is a histogram of A<J>. The standard deviation, 0 , , of the 

A<p 

distribution is 36°. This variation reduces measured field-aligned 
anisotropies. An estimate of the size of the reduction is obtained by 
approximating the distribution of A<J> as a Gaussian distribution with 
a standard deviation of 36°. Such a distribution of field directions 
reduces the measurement of field-aligned anisotropies by a factor of 
0. 82 , 


The 6-hour average magnetic field directions are later com- 
pared to the diffusive anisotropy direction. The variation in the 
magnetic field direction is a sc \le size for the root mean square dif- 
ference between the computed field direction and the effective field 
direction seen by the average particle. For example, occasionally 
only 5 hourly averages are available to compute the 6-hour average 
field direction. The root mean square difference between the com- 
puted field direction and the direction that would have been com- 
puted had all 6 hourly averages been available is ~ a , / 6 or — - 6° . 

A <p 
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Figure IV-10 

The lower panel is the distribution of the differences in the 
hourly average magnetic field direction * and the 6-hour average 
field direction ? in which the hourly average is included. 

The upper panel is a histogram of the number of hourly averages 

included in the 212 6-hour periods. 
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V. DATA SELECTION 

A . Intr oduction 

The primary data set for this study runs from 72/273 to 74/2 . 
The set begins when the first useable information from the IMP-7 
satellite was obtained; the end of the set is the last time for which 
simultaneous solar wind velocity and magnetic field direction have 
been obtained. Figure V-l is a plot of the PLO rate for this time 
span. Those 6-hour periods used for the results of this study having 
a plO rate a 0. 01/second are identified. Some periods are not used 
because either the proton or plasma data are not available, while 
others have been rejected in order to avoid possible sources of error. 
In addition, since the purpos * this study is to investigate stream- 
ing during periods between prompt events, periods dominated by 
prompt events are not included in the primary data set. For much of 
the analysis a minimum proton flux is required in order to avoid pe- 
riods when computed anisotropies have large statistical uncertainties. 
An analysis of the selection process and of the possible sources of 
error that have been avoided follows. 

As indicated in the previous chapter, 6 -hour averages of the 
proton and plasma data are used. There are 1839 6-hour periods in 
the time span from 18:00 UT 72/273 to 12:00 UT 74/2. Proton fluxes 
and their associated anisotropies are available for 1720 (94 per cent) 
of these periods, magnetic field data for 1543 (84 per cent) periods, 
and solar wind data for 1586 (86 per cent) periods. There are 1264 
(69 per cent) periods for which all these data are available. 
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Figure V-1 

The 6-hour averages of the PLO rate for the time span used in this 
work. Periods when the satellite is sunward of the earth are indi- 
cated. Periods when the PLO rate is between 0.01 /second and 0.3/second 
that are used in this work are also indicated. 
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B. Instrumental Effects 

Several possible instrumental sources of error were discussed 
in Chapter IY. Bong term biases affecting the anisotropy of sectored 
rates were shown to be small by computing the anisotropy of two 
rates dominated by an isotropic background. It was shown that 
anisotropies using sectored rates and analyzed events are consistent. 
An upper limit was placed on the background rate for D2 events which 
indicated that the anisotropy of D2 events is negligibly affected by 
any such background. The effect of these possible sources of errors 
is sufficiently small so that no periods have been eliminated to avoid 
them. Very rarely an error is made in the sectored rate data that 
produces a very large number of counts in one sector. These periods 
have been identified by the poor fit made to the data when computing 
the anisotropy. There are nine such 6-hour periods, each having a 
goodness -of-fit parameter \ v greater than 200. There is no indica- 
tion of any similar effect in the analyzed events or other sectored 
rates during these periods. Nevertheless, no data from these 9 pe- 
riods when this malfunction occurred are used in this study. 

Large fluxes can introduce errors into anisotropy measure- 
ments of sampled events such as D2 analyzed events. A bias is put 
into the anisotropy measurement because the instrument can transmit 
only one analyzed event every 0. 64 seconds. Sectored rates are not 
affected by this problem and so can be used to estimate the size of 
errors. In the limit of large fluxes the analyzed events become 
isotropic. Roelof (1974) has made a theoretical study of the effect 
for a class of instruments similar to the Caltech experiment. The 
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major difference in the Caltech experiment is that the priority assigned 
to events during the current readout period depends on whether the 
most recent event read out was a narrow geometry electron-type event. 
However, this changing priority structure has a small effect on the 
measured anisotropy of D2 events during the periods used in this 
study. The previous readout is usually a narrow-geometry proton 
(rate ~ PLO), a wide-geometry event (rate >0.3 second), or a neutral 
event (rate ~ 0. 1 /second) rather than a narrow-geometry electron- 
type event (rate * ELO + EHI ^ 0. 1/second). So to a first approxi- 
mation, narrow- geometry electron events always have the highest 
priority. This reduces the live time for D2 analyzed events by the 
probability of getting a narrow-geometry electron-type event during a 

readout period 6 per cent. An anisotropy in the electron rate 

would make the live time depend on the spacecraft viewing angle and 
thus affect the anisotropy measured using D2 analyzed events. As 
discussed in Chapter IV, ELO and EHI are usually isotropic. The 
anisotropy for ELO + EHI rates summed over all periods used in this 
study with 0. 01/second < PLO < 0. 03/second is 0. 5 per cent ±0.4 
per cent. Thus, the difference in live times introduces an error of 
~ (°- 5 P er cent) * (6 per cent), or 0.02 per cent. 

Since the history dependence of the priority system of the Cal- 
tech instrument has a minimal effect, the rate dependence of the an- 
isotropy of D2 analyzed events is compared to the formulas derived 

by Roelof for small anisotropies: 

1^1 i 

r 3 TT^T = — r [( 1+n2 — os ^ - 1 HnVl-nliS. Y . +n - ncos Y \f 

(5-1) 
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tan(A<|>) = ta n(<i> OBS -‘! , trU e ) " n 


e n ^ - i -n(siny - n cosy-fn) 
e 11 ^ -1+n (n siny+cosy - 1 ) 


n = mean event rate/ radian 


(5-2) 


Y = (To-TjJW 

U) = rotation rate = 4.77 radians /sec 

t = read out period = 0. 64 seconds 
o 

= dead time for read out = 0. 03 seconds 
The anisotropy of the PLO rate is used as a measure of the true or un- 
biased anisotropy of D2 analyzed events. As noted in Chapter IV, 

PLO events have a slightly larger energy interval than D2 analyzed 
events, and so the anisotropies of the two event types may be only ap- 
proximately equal. 

The mean value of r and A4> are plotted as a function of the 
PLO rate in Figures V-2 and V-3. The solid lines are the values 
predicted by eqs. (5-1) and (5-2) in the approximation that all priority 
0 and 1 events are PLO events. Only periods when both anisotropies 
are statistically significant ( £/<^ >2.5) are included. 

The observed values of r and A4> are consistent with the pre- 
dicted values except r does not approach 1 for small PLO rates. 
However, r is essentially constant for PLO less than 0.3/second, 
indicating that the discrepancy is not a biasing effect caused by hign 
rates. Consequently, only periods when the PLO rate is less than 
0. 3/second are included in the final data set of this study. At a PLO 
rate of 0. 3/second, eqs. (5-1) and (5-2) give r = 0. 95 and A$ = 3 
with correspondingly smaller values for lower rates. Thus, any 
biasing effect caused by high rates will minimally affect the observed 
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Figure 7-2 

The mean of the ratio of the anisotropy amplitude determined using 
D2 analyzed events to the anisotropy amplitude determined using 
the sectored PLO rate as a function of the PLO rate. Only periods 
when the anisotropy amplitude e is more than 2.5 times the statistical 
uncertainty in the anisotropy amplitude for both D2 and PLO 
anisotropies are included. Periods when the average ELO rate is 
larger than 0.15/second are not used. The ratio calculated using 
eq. (5-1) is shown by the solid curve. 


anisotropy amplitude 

PATE DEPENDENCE 


72 / 273 - 75/15 
-90° < <£ sat < 90 
ELO< 0.15 



PLO (counts/second) 
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Figure V-3 

The mean difference in the direction of the anisotropy determined 
using D2 analyzed events and the direction of the anisotropy 
determined using the sectored PLO rate as a function of the PLO 
rate. The difference calculated using eq. (5-2) is shown by the 


solid curve. 


ANISOTROPY DIRECTION 
RATE DEPENDENCE 


PLO (counts/second) 
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anisotropies for the periods used in this study. 

Since § = 3S/wU , an estimate of the differences of anisotro- 
pies using FLO and D2 analyzed events is given by the ratio of the 
average particle velocities. The ratio of 0. 96 for protons is not 
quite as small as the observed ratio of - 0. 91 at low rates. The re- 
maining factor of ~ 0. 95 is not understood, but is a small effect. 
Typical anisotropies used in this study are - 10 per cent. A correc- 
tion factor of 0. 95 would mean the true anisotropy is 10. 5 per cent. 
Because the factor of - 0. 95 is not understood, the anisotropy of D2 
analyzed events is used as measured. 

No attempt has been made to correct the Observed anisotropy 
of D2 analyzed events for high rates; instead, periods when the PLO 
race is larger than 0. 3/second are not included in this study. Equa- 
tion (5-1) indicates that the measured anisotropy will be 0. 95 of the 
true anisotropy at a FLO rate of 0. 3/second. This typically is an 
absolute error of — 0. 5 per cent In the anisotropy measurement. 

Most periods used have smaller rates and so correspondingly smaller 
errors . 

Anisotropies could be determined at higher fluxes using FLO 
sectored rates. This has not been done since periods with FLO rates 
> 0. 3/second are predominantly due to prompt solar events, while 
the purpose of this work is to study the periods between prompt 
events. 

C. Magnetospheric Influence 

The purpose of this work is to study the interplanetary stream- 
ing of protons. Consequently, only periods when the satellite is out- 
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side the earth's bow shock -- the first interaction of the earth's mag- 
netic field with the solar wind -- are included in the final data set of 
this work. Data taken when the satellite is inside the bow shock are 
used only in this section for comparison with data taken in inter- 
planetary space. 

To limit the data set to interplanetary space, only periods 
when the satellite is sunward of the earth are included. This corre- 
sponds to a range from -90° to +90° in the Geocentric Solar Ecliptic 
longitude of the satellite position. As indicated in Figure II- 1, this 
is outside of the average position of the earth's bow shock as deter- 
mined by Behannon (1968) during 1966. 

The presence of the earth could influence observed anisotro- 
pies even when the satellite is outside the bow shock. The magneto- 
sphere could be a source of low-energy protons or could distort the 
local flow of interplanetary protons. Krimigis, et al. (1975) have re- 
ported a highly anisotropic flow of 0. 29 to 0. 5 MeV protons from the 
direction of the earth using data from the JHU/AFL experiment a- 
board IMP-7 while IMP-7 was in interplanetary space. Using the 
Caltech experiment on IMP-7, Mewaldt, et al. (1975b), looking at 
similar periods , did not find this streaming at the higher energy in- 
terval of 1. 3 - 2. 3 MeV -- the same energy interval used in the pres- 
ent study. The absence of streaming for particles > 1 MeV is con- 
sistent with the source of these particles proposed by Hovestadt and 
Scholer (197 6) to explain the source of the lower energy fluxes. 

The size of any possible magnetospheric effect on interplane- 
tary particle streaming is determined by the interaction of the mag- 
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netosphere with the interplanetary medium as seen by ~ l| MeV pro- 
tons, which have a Larmor radius of ~ 5 earth radii. Due to the 
complexity and variability of this interaction, no attempt has been 
made to calculate theoretically the effect of the magnetosphere. 
Rather, the dependence of the observed anisotropy of the satellite po- 
sition has been examined, A dependence of the anisotropy on satel- 
lite position would be expected if the earth blocked the flow of par- 
ticles from beyond 1 AU. This blockage would be more effective 
when the satellite is near -45° when the field line connecting the 
satellite to regions beyond 1 AU must drape over the magnetosphere 
than when the satellite is near +45 when the corresponding field line 
usually will not encounter the magnetosphere. 

Figure V -4 shows the anisotropy amplitude observed as a 
function of spacecraft position. To reduce statistical uncertainties 
in the individual points, 6 -hour averages have been combined into 
daily averages. An outstanding feature of the plot is the compara- 
tively small anisotropies seen when the satellite is near 180° -- 
roughly when the satellite is inside the magnetopause. Outside of ~ 
140° to 210° the anisotropies are at least qualitatively independent of 
satellite position. For more quantitative comparisons, the data have 
been grouped into 3 regions: 

Region 1 0° to 90° interplanetary space 

Region 2 -90 to 0° interplanetary space 

Region 3 140 to 210 — > magnetotail 

The final data set includes data from only Regions 1 and 2. If the 
magnetosphere has little influence on anisotropies in interplanetary 
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Figure V -4 

The amplitude of the observed anisotropy as a function of the longi- 
tude of the IMP -7 satellite in Geocentric Solar Ecliptic coordinates. 
A typical ± Iff error bar is indicated. 
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72/273-75/18 
0.03 < PLO < 0.3 
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space, Region 1 and Region 2 will have similar distributions. Histo- 
grams of the observed anisotropy amplitudes are shown in Figure V- 
5 for each of the above regions. The mean anisotropies for the re- 
gions are; 

Region 1 12. 6 per cent ±1.7 per cent 

Region 2 13. 5 per cent ± 1. 0 per cent 

Region 3 3.8 per cent ± 0. 5 per cent 

The difference of the mean anisotropy of Regions 1 and 2 is not sta- 
tistically significant; the Student's t-test indicates more than a 60 
per cent probability of observing as large a difference due to chance 
alone. In contrast, the mean anisotropy for Region 3 is only 1. 2 
standard deviations above the mean amplitude Gf 3.2 per cent expect- 
ed if the flux were isotropic. Thus, while the magnetosphere has a 
significant effect on anisotropy measurements when the satellite is in 
the magnetotail, there is no indication of magnetospheric influence 
on measurements taken in interplanetary space. 

The magnetosphere could affect the observed anisotropy but 
still maintain the same average amplitude. Figures V-6 and V-7 
plot the individual observed anisotropy vectors for Regions 1 and 2, 
respectively, using 6 -hour averages. Again, the distributions for the 
two regions are similar. The characteristics of the distributions of 
the x- and y- components are: 

Region 1 Region 2 


No. 

Mean 

Sigma 

Median 


151 
-7.82 # 
13.4 # 

-8. 12 $ 


-2. 15 £ 
14. 6 $ 

-3.49 $> 


205 

-9. 12 $ 
11.9 $ 

-9. 39 $ 


-3. 90 $ 

12 . 0 <f> 

-4.63 # 


Figure V-5 

Histograms of the amplitude of the observed anisotropy for three inter 
vals of the longitude of the IMP-7 satellite. The mean amplitude is 
indicated for each histogram. 
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Figure V-6 

The observed anisotropy for periods when the longitude of IMP-7 is 
between -90° and 0°. The data are from 72/273 to 75/18. The means of 
the x- and y-components of the anisotropy are indicated. 



Figure V- 7 

The observed anisotropy for periods when the longitude of IMP-7 is 
between 0° and +90°. The data are from 72/273 to 75/18. The means 
of the x- and v-components of the anisotropy are indicated. 
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The result of applying the t- test to the differences in the means is 
that there is no significant difference in either component average at 
the 80 per cent confidence level. The difference in the sigmas can be 
compared using the F-test. The probability of finding differences as 
large as those observed due to random fluctuations is 5 per cent for 
the x-component and 1 per cent for the y-component. 

The apparently significant difference in the spread of the dis- 
tributions must be interpreted with caution. The difference reflects 
the existence of a few periods with anisotropies much larger than 
typically found. For example, removing two periods from Region 1 
will make the sigma in the x-direction smaller than that found in Re- 
gion 2. Thus, the estimate of the chance likelihood of such a differ- 
ence using the F-test appears too small. This reflects the existence 
of non-Gaussian tails of the observed distributions. 

In summary, the evidence indicates that any magnetospheric 
influence on the anisotropies observed in interplanetary space is 
small. No significant difference was found in either the mean ampli- 
tude, the mean x-component, or the mean y-component of the aniso- 
tropy when measured on both sides of the sun - earth line. A differ- 
ence of about 2 per cent was found in the spread of observed anisotropy 
components, but this could be due to the chance occurrence of two 
periods of large anisotropies. Consequently, data from both Region 1 
and Region 2 are included in the final data set. 

D- Prompt Solar Events 

In order to study periods between prompt events, prompt events 
must be identified and eliminated from the data set. This is done in 
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two ways. First, periods when the PLO rate is more than 0. 3/sec- 
ond (a proton flax of 1. 2/cm -sec-sr-MeV) are not included. This 
eliminates periods dominated by large prompt events. Second, indi- 
\ idual prompt events are identified and days near the beginning of the 
observed increase are eliminated. Prompt events are identified by 
their sharp increases in flux, typically a factor of 10 in 12 hours. 
Corroborative characteristics are large peak fluxes, velocity dis- 
persion, exponential decay, and an accompanying prompt electron 
event. Anisotropy data are not used. Table V-l lists the 25 prompt 
events identified for the period 72/27 3 to 74/2. The day preceding 
the sharp increase, the day of the sharp increase, and the two days 
following the sharp increase are not used in this study. Only whole 
days are eliminated. Rejecting these periods eliminates the begin- 
nings of prompt events e\ en though the PLO rate may not have 
reached 0. 3 /second. Small prompt events are also eliminated. The 
day previous to the sharp increase is eliminated to avoid possible 
srnall injections preceding the main injection of particles at the sun. 


YEAR 

1972 

1973 


i 
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Onsets of Prompt Events, 72/273 - 74/2. 

DAYS 

282, 291, 303, 329, 333, 348, 351, 

46, 71, 78, 89,102,114,119,154,171,180,210,250, 
261, 270, 277, 292, 307, 310, 
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VI. RESULTS 

A. Introduction 

The observed anisotropy of D2 analyzed events has been de- 
termined according to the procedure discussed in Chapter IV. D2 
analyzed events are nominally 1.3 to 2.3 MeV protons. The periods 
used have been selected according to the criteria discussed in Chapter 
V. The criteria eliminate periods dominated by prompt solar events. 

Both observed and diffusive anisotropies are used. The diffu- 
sive anisotropy was defined in Chapter III: 

t P 3CV 

*DIF ^OBS " w 

For small anisotropies, the diffusive anisotropy is the anisotropy in 
the rest frame of the solar wind (Forman, 1970; Balogh, et al. , 1973), 

The observed anisotropy is displayed using the Solar Ecliptic 
coordinate system. The diffusive anisotropy is displayed using both 
the Solar Ecliptic coordinate system and what will be termed the Mag- 
netic Ecliptic system. The Magnetic Ecliptic system is obtained by 
rotating the Solar Ecliptic system about its z-axis until the new x-axis 
is colinear with the direction of the projection of the observed mag- 
netic field onto the ecliptic plane. The Magnetic Ecliptic system is 
used to display the importance of the magnetic field direction on the 
diffusive anisotropy. Because particle propagation is independent of 
the sense of the magnetic field, there is an arbitrary choice in the 
sense of the x-axis. The positive x-axis is defined to be within 90° of 
315° (the long-term average magnetic field direction) in Solar Ecliptic 
coordinates. The y-axis is chosen to form a right-handed coordi- 
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nate system. The new x- and y- axes are referred to as *, and 

K± reSPeCtiV6lyi the Pr ° jecti0n of the diffusive anisotropy onto these 
axes are designated and g A (see Figure VI- 1). 

The Magnetic Ecliptic coordinate system is constructed so 
that particles flowing along the magnetic field lines toward the sun 
nominally produce a positive S| . However, when the observed field 
direction is nearly perpendicular to the long-term average field di- 
rection, it is not certain which direction along the field leads to the 
sun. The direction nearest is used, but may not be correct. If in- 
correct, a positive U would correspond to flow away from the sun. 

Of the periods used when PLO is between 0. 01/second and 0. 3/sec- 
ond, only 7 percent have field directions differing by more than 75° 
from the long term average. 

^ — 6^Hour A verage Anisotropy 

The observed anisotropies are plotted in Figure VI-2. The 
diffusive anisotropies for the same 112 periods are shown in Figures 
VI-3 and VI-4 using the Solar Ecliptic and Magnetic Ecliptic coordi- 
nate systems respectively. The PLO rate for each period included is 
greater than 0. 03/second and less than 0. 3/second. 

The mean values for each of the groups are: 


(i) 


<*x> 

v 


b OBS 
S. E, y 

-7. 0*0. 5 2.4 

-1.7*0. 5 2.2 


^DIF ^ 


S.E. y 

13.6*0.5 2.5 

-3. 1*0. 5 2.2 


M.E. x 

v 

10.7*0.5 3.2 

1.3*0. 5 2.0 


Weighted means are used - each data point is weighted according to 
the statistical uncertainty in the anisotropy measurement. The error 


Figure VI-1 


The upper panel shows the ecliptic plane as viewed from the North 
Ecliptic Pole. The x-axis and the y-axis of the Solar Ecliptic coor- 
dinate system are shown. The bottom panel shows the rotation of the 
Solar Ecliptic coordinate system about the z-axis used to define the 
Magnetic Ecliptic coordinate system. The axes of the Magnetic Eclip- 
tic system are parallel to and perpendicular to the projection of the 
magnetic field onto the ecliptic plane. 
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COORDINATE SYSTEMS 



Figure VI- 2_ 

The observed anisotropy, f 0BS > using Solar Ecliptic coordinates. Each 
dot represents one 6-hour period. A typical ±lo error bar is indicated 
as well as the means of the x- and y-components of the anisotropy. 



Figure VI -3 


The diffusive anisotropy, l DI p» using Solar Ecliptic coordinates. Each 
dot represents one 6-hour period. A typical ±1it error bar is Indicated 
as well as the means of the x- and y-components of the anisotropy. 




Figure VI-4 

The diffusive anisotropy, using Magnetic Ecliptic coordinates. 

Each dot represents one 6 -hour period. A typical ±lc error bar is 
indicated as well as the means of £u and . 
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»™onnnntBlLirY OF THE 


bar 8 quoted on the means are determined from the statistical uncer- 
tainties in the individual measurements. The \^b being larger than 
1. 0 indicates that the typical statistical uncertainty in each measure- 
ment of ~ 5 per cent is smaller than the real variation in the observa- 
tions. Thus, a better estimate of the uncertainty in the mean is given 
by the product of the quoted error bars and • The convention is 
adopted that unless i he °f the fit is explicitly stated, quoted error 

bars will have been multiplied by . 

The typical SqbS - 8 awa y from the sun in a direction slightly 
counterclockwise from radial. The typical radial component is 7 
per cent. The diffusive anisotropy is qualitatively different from the 
observed anisotropy with the typical flow back toward the sun with a 
radial component of ~ 14 per cent. The mean £|j is somewhat small- 
er than the mean 6 DIF x - This is due to the larger number of peri- 
ods having a negative £u than having a negative x • 

Equation (3-4) relates the diffusive anisotropy to gradients in 
the particle density. The typical diffusive flow back toward the sun 
indicates that the typical radial gradient is positive --a larger den- 
sity beyond 1 AU than inside 1AU. In fact, 86 of the 112 periods (77 
per cent) included in Figure VI-4 have a positive £ || . Two effects 
could systematically bias this percentage. First, large statistical 
uncertainties in the individual measurements tend to make the periods 
evenly split between positive and negative £j| . However, if the dis- 
tribution of values around the mean of 10.7 per cent were all due to 
the typically 5 per cent statistical uncertainty in the individual meas - 
urements, only 2 of the periods would have a negative §jj . The other 


possible effect is the mis -identification of which direction along the 
magnetic field leads toward the sun. Such a mis -identification as- 
signs the wrong sign to £|). This mis -identification is most likely 
when the field is far from its long-term average direction. A larger 
fraction of the periods with negative £ n than those with positive £|j 
have magnetic field directions more than 45° from the long-term av- 
erage direction. Sixteen of the 26 periods with flow away from the 
sun have such a field, and 16 of the 86 periods with flow toward the 
sun have such a field. Of the 32 periods with a magnetic field direc- 
tion more than 45° from the long-term average direction, 12 either 
are associated with magnetic field sector boundaries or have field di- 
rections opposite the current sector’s prevailing field direction. 
Eliminating these 12 periods from consideration increases the 77 per 
cent quoted above to 82 per cent. A nominal value of 80 per cent will 
be used in Chapter Vll. Eliminating these 12 periods also increases 
(£ll> and <£ x ) to 12. 1 per cent and 1. 8 per cent, respectively. These 
increases are comparable to the statistical uncertainties in these 
values. 

C, The Diffusive Anisotropy 

The dependence of the diffusive anisotropy on various parame- 
ters is now investigated. From eq. (3-4), 

fbi F = <-!>»' (VU/TX) . (6-1) 

Although & and VU are not known, eq. (6-1) does suggest what quan- 

— V 

tities might affect Spgp • 


1. Flux Dependence. Figure VI- 5 contains histograms of § j| 
for 4 groups of data selected by the PLO rate. PLO rates from 
0. 00 3 /second to 0* 3 second are included. The weighted means of the 
distributions are: 



PLO 

periods 

<g||> (*) 

±(t) 

0 . 1 

-0.3 

24 

5.4 

2.8 

0.03 

- 0 . 1 

88 

11.9 

1.9 

0.01 

- 0.03 

100 

10.0 

1.7 

0, 003 

- 0.01 

105 

15.4 

2.6 


The weighted mean for the 4 groups is 10. 8 per cent. There is 
a 5 to 10 per cent chance of the variations from the mean being as 
large as observed due to statistical fluctuations. Thus, there may be 
a statistically significant flux dependence to the anisotropy. However, 
the typical diffusive streaming is back toward the sun for all the 
groups, indicating streaming toward the sun to be typical for a wide 
range of low-energy proton intensity. 

2, Time Since Most Recent Prompt Event . Prompt solar 
events can dominate low-energy proton fluxes at 1 AU for several 
days. Periods dominated by prompt events have been eliminated from 
this study by the criteria discussed in Chapter V. Prompt events 
could have some influence on gradients and thus diffusive anisotropies 
for longer periods. For example, McKibben (1972) has reported de- 
cay times as long as 87 hours late in some prompt events after an 
initial decay time of ~ 24 hours. 

To investigate possible influence, ||j is plotted in Figure VI-6 
as a function of the elapsed time, At, since the onset of the most re- 
cent identified prompt event. The onset times of prompt events have 



Figure VI-B 

The distributions of g, for four intervals of the PLO rate. The typical 
±lo error bar for an individual measurement Is indicated for each group. 
The mean g n for each group is also indicated. 
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Figure VI-6 

The dependence of C« on the number of days since the onset of the most 
recent Identified prompt event. A typical ±1* error bar Is shown. The 
dashed line Is the least squares fit to the data. 


I 


TIME SINCE PREVIOUS PROMPT EVENT 

72/273-74/2 
0.03 < PL0<0.3 
— 90° < <f> < 90° 


# e„ = (6.l±0.9)+ 
(0.30±0.05)At 


At (DAYS) 
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been listed in Table V-l. The scatter in the data appears larger than 
any trend, indicating that At is not a dominant factor in determining 
£ |i . A linear fit to the data: 

|||($) = (6. 1 ±0.9) + (0. 30 ± 0. 05)At (6-2) 

X = 3.2 
v 

where At is in days, shows a small trend toward a larger anisotropy 

for larger At, although the large x v means the fit is poor. The line 

determined by eq, (6-2) is positive for all values of At greater than 

0 , indicating that the typical diffusive anisotropy is toward the sun 

for all At included in the data set. 

3. Magnetic Field Direction. As discussed in Chapter III, 

the magnetic field is understood to control the diffusion of particles in 

interplanetary space and thus affect f DJF . For example, if (kj./ku) 

« 1 , then £ dif will be parallel to the magnetic field direction, <k_. 

B 

On the other hand, if Hi = H|| (isotropic diffusion), then should 

DIF 

be independent of . 

To investigate the relative size of h± and k \\ , the direction of 

the diffusive anisotropy, 4> mF , has been plotted as a function of <J> B 

and a linear fit made to the data. Both the field and the anisotropy 

are 6-hour averages. The fit is complicated by the independence of 

the streaming on the sense of the field direction. Consequently, <> 

T B 

has been restricted to a 180 range by the following mapping: 



*OBS 


for 

- 90 ° < *OBS < 90 ° 

II 

£ 

*OBS + 

180° 

for 

*OBS < - 90 ° 


*OBS " 

180° 

for 

*OBS > 90 ° 


where 4>Qgg is the observed field direction. 
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A linear relationship between 4> DI p 

W = a + b 


and 4> is assumed: 

D 

= f(4> B ) . 


If diffusion is isotropic, b is expected to be consistent with 0. For 

field -aligned , a = 0 $ ad b = 1 . The parameters of the fit -- 

XjXjT 

a and b - - are determined by minimizing 

x 2 = r r DIS ^DiF.i> _ 

^DIF, i 

where i indexes the individual data points and the distance measure is 
given by 


DET(* 4^) = (4> r ^) 

mod 180 


(6-3) 


and the range of DIST is [-90°, +90°). The statistical errors in the 

2 

parameters are estimated by the change required to increase x by 1 . 
The optimal parameters were determined by a grid search technique. 
The best fit is 


<j> DIF = (-10. 2° ±1.3°) + (0.71 ±0.03)4.3 (6-4) 

with x =2.2 

Figure VI-7 contains two panels. One panel is of <W versus 4> B i 
in the other panel, the diffusive anisotropy directions have been re- 
plotted to show their distances [eq. (6-3)] from the best fit [eq. (6-4)]. 

The strong dependence of the diffusive anisotropy direction on 
■die magnetic field direction shown by eq. (6-4) indicates that diffusion 
is not isotropic, so that kj. is less than Ky - The slope of the line de- 
termined by eq. (6-4) being less than 1 suggests that kj. is not negli- 
gibly small. The value of (x.i/H|| ) consistent with the data depends on 
the relative gradients parallel to and perpendicular to the field lines. 


0 


Figure VI-7 

The left panel shows the direction of the diffusive anisotropy as a 
function of the magnetic field direction. The data are replotted in 
the right panel and the least squares fit Is indicated by the dashed 
line. Points which would have been plotted outside the solid lines 
have had +180° or -180° added to their anisotropy direction. A typical 
±la error bar in the anisotropy direction is Indicated. 


72 / 
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The determination of the value will require further study. 

4. Magnetic Field Fluctuations . As seen in eq. (6-1), the 
size of the diffusion tensor could affect • T ^ e work of Jokipii 

(1966) relates x to fluctuations in the magnetic field with wavelengths 
comparable to the particle's Larmor radius. As discussed in Chapter 
III, the exact relationship is still a subject of controversy. This 
study has sufficient information only to estimate the amount of tur- 
bulence in the magnetic field over a wide range of frequencies. The 
magnetic field data used in this work are based on hourly averages 
provided by the National Space Science Data Center. These hourly 
averages are in turn derived from averages of even finer time scales 
which range from 1. 3 to 48 seconds, depending on which experiment- 
er's data are used. The standard deviation in the fine time scale 
averages used to determine the hourly average field is a measure of 
the fluctuations in the magnetic field with frequencies between 
~ (1 sec)” 1 to -(1 hour)” 1 . The frequency appropriate for l| MeV 
protons is ~ (70 sec)" 1 . A dimensionless variable is given by B/tfg 
for each hour. Figure VI-8 is a plot of the average of the hourly 

B/ct versus the amplitude of the diffusive anisotropy, and the best 
B 

linear fit to the data. In view of the large , the difference of the 
slope from 0 is not considered significant. 

This result is consistent with a diffusive anisotropy independ- 
ent of the magnitude * , but the necessarily crude method for esti- 
mating H prevents any definitive determination. 

5 , S olar Wind Velocity. The solar wind velocity is an im- 
portant parameter in particle propagation, so changes in the velocity 


Figure VI -8 

The dependence of the amplitude of the diffusive anisotropy on the 
relative fluctuations in the magnetic field. A typical ±la error bar 
in the amplitude of the anisotropy is indicated. The dashed line is 
the least squares fit to the data. 
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B FLUCTUATION 


72/273-74/2 
0.03< PL0<0.3 
-90° < <f> sat <90° 
NO PROMPT EVENTS 


II 






lei = (21 ± i) 

+ (- 1.0± 0.3)(B/Oo) 
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V 


1 


could alter spatial gradients and thus the diffusive anisotropy. Figure 
VI -9 plots for 4 groups of solar wind velocity. The components 


are: 



V gw (km/sec) 

Periods 

<! x > <*) 

< 400 

87 

7.8 ± 1.5 

400 - 500 

69 

15.7 ±1.5 

500 - 600 

48 

14.7 ±1.6 

> 600 

8 

17.8 ± 3.8 


The extreme velocities used are 270 km/sec and 660 km/sec. 
The variation in the mean is statistically significant for the x -compo- 
nent, indicating an increased streaming toward the sun for periods 
with large solar wind velocities. A model consistent with such a cor- 
relation will be presented in Chapter VII. 

D. Particle Flow 

As noted in Section VI-A,the diffusive anisotropy is the aniso- 
tropy in the rest frame of the solar wind for small anisotropies. The 
particles are also being convected with the solar wind, so the aniso- 
tropy determined by the actual particle flow is given by 

-* -* 3 y 

^part = ^DIF + w~ ' 

This differs from the observed anisotropy in that the convective term 
is not multiplied by the Compton - Getting factor. The Compton - 
Getting factor corrects for the difference in energy as a function of 
direction and so is not appropriate for the actual flow of particles. 

Figure VI-10 is a histogram of the x-component of £p art * 

The mean value, weighted by the intensity PLO, is 


Figure VI-9 

The diffusive anisotropy, ? D jp» using Solar Ecliptic coordinates for 
periods grouped by the solar wind speed. Typical ±l<r error bars are 
shown. The means of the x- and y-components of the anisotropy are 


indicated 








72/273-74/2 
0.01 <PLO <0.3 
500 < V sw < 600 

-90° < </>sat <90° 



72/273-74/2 
0.01 <PLO< 0.3 
600 < Vgyy < 700 

-90° < <£sAT < 90° 
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Figure VI -10 

The distribution of the x-component of t part j cle * Positive values 
indicate particle flow toward the sun. The mean value, determined 
by weighting the individual values by PLO, is indicated. 
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PARTICLE ANISOTROPY HISTOGRAM 


p 72/273-74/2 
| O.OI < PL0<0.3 
i -90° < </> SAT < 90° 



) 0 50 
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U = 


S(£ . .«PLO) 

P ar *’ * = 5.0$ ± 0.5$ . 


D PLO 

The average particle streaming is 

4tt <K„ a ,t v ) pLO > 2 

(S ) = P * — - — _ — r - .052 ±.005 particles/cm -sec 

N part, x 3 (geometric factor) r 

between 1. 2 and 2. 4 MeV. 

At an average particle energy of 1. 6 MeV, this corresponds to 

2 

an energy flux of 0. 083 MeV/cm -sec flowing toward the sun. 

For comparison, the energy flux in the bulk flow of the solar 

c 2 

wind is ~ 2 X 10 MeV/cm -sec. The energy flux in the interplane- 
tary magnetic field is ~ 1 per cent of this value. Thus, the energy 

_7 

flow of 1.2 to 2. 4 MeV protons at 1 AU is only — 4 X 10 times as 

-5 

large as the solar wind bulk flow energy and ~ 4 X 10 times as 
large as the energy flux in the interplanetary magnetic field. 

Interplanetary acceleration, perhaps at the interface of high- 
and low -speed solar wind streams, has been suggested by McDonald, 
et al. (1975 ) to explain increases in low-energy protons seen near 3 
AU. Such acceleration could also be the source of the particles ob- 
served in this study. The relatively small energy flux of ~ 1.6 MeV 
protons seen at 1 AU is consistent with the solar wind being the 
energy reservoir for the acceleration. A definitive energy balance 
study must await knowledge of the streaming at radial distances be- 
yond the source of the particles as well as knowledge of the specific 
acceleration mechanism. 


118 


VII. DISCUSSION 


REPRODUCIBILITY of th 

. _ -^.1 n-n TC? DA( r. 


A. Source 

1. Spatial Gradient. The Caltech EIS experiment aboard 
IMP-7 has been used to measure the streaming of cosmic -ray protons 
in the energy interval 1. 3 -2. 3 MeV during selected periods from 
September 1972 to January 1974. As discussed in Chapter V, the 6- 
hour periods were selected to study the times between prompt solar 
particle events. Statistically significant anisotropies could not be 
measured for individual periods at the lowest flux levels observed, so 
the periods used are during enhancements of the flux above quiet 
time levels. For 112 6-hour periods when the PLO rate was between 
0. 03 and 0. 3/second, the amplitude of the average was ~ 7 per 

cent. The particles were observed to be streaming from the sun in a 
cirection slightly counterclockwise from radial. The diffusive aniso- 
tropy has been determined by subtracting the convective anisotropy 
from the observed anisotropy. The average ? DIF Bad an amplitude 
of ~ 12 per cent, with the flow toward the sun along the observed 
magnetic field direction. The flow along the field was toward the sun 
in ~ 80 per cent of these periods. This flow back toward the sun was 
found to be typical of all the PLO rates used in this study - - from 
0. 003/second to 0. 3/second. 

The previously proposed models of delayed events seen near 
1 AU, which were reviewed in Chapter I, involve continuous injection 
of ~ 1 MeV particles into the interplanetary medium near the sun, 
thereby producing flow away from the sun. 

Flow back toward the sun implies a positive radial gradient in 


T ‘ T r 


i'- “,V 1 ., 1 .W* » - 1 11 
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the particle density. Evidence of the existence of such a radial gradi- 
ent, at least during certain periods, has been presented by McDonald, 
et al. (1975 ) on the basis of the average relative size of delayed 
events seen at 1 AU and at ~ 3 AU from November 1973 to April 1974. 
Whether the delayed events seen near 3 AU and the increases of the 
present study are the same phenomena has yet to be determined. The 
determination is made difficult by the small overlap in observing time 
and the separation in azimuthal position of the two spacecraft. In any 
case, some of the models proposed to explain the increases near 3 AU 
may be relevant to the anisotropies reported herein. 

2. Interplanetary Acceleration . McDonald, et al. (1975) 
considered interplanetary acceleration the most likely explanation of 
the growth of delayed events with increasing radial distance. Fisk 
( 1976 ) has developed numerical solutions to the Fokker -Planck 
propagation equation involving interplanetary acceleration to explain 
the observed increases. The particles are injected at low energies 
( « 1 MeV) near the sun and are accelerated as they propagate 
through the interplanetary medium. Acceleration is assumed to take 
place throughout the solar cavity. The rate of acceleration is ad- 
justed to produce an increase in the intensity by a factor of ~ 10 be- 
tween 1 and 3 AU. Two acceleration processes are considered. The 
rate of the first process, Fermi-acceleration, is determined by the 
diffusion coefficient k|| ; a H|| of - 1. 8 X 10 19 cm 2 /sec was needed 
to produce the desired radial gradient. Such a H|| produces a diffu- 
sive anisotropy at 1 AU of ~ 0. 25 per cent. This value is ~ 50 times 
smaller than the typical anisotropy found in the present work. Thus, 



the model is not consistent with the observations. The second accel- 
eration process, transit-time damping, does not depend on K|j and so 
cannot be tested by the anisotropy results reported in this work. 

Acceleration may be taking place in localized regions rather 
than throughout the solar cavity. Pesses, et al. (1976 ) have pre- 
sented evidence that the acceleration of protons in co -rotating shocks 
at the boundary between two solar wind streams is a common feature 
of interplanetary space at distances ^2.6 AU. Propagation along the 
boundary could spread the azimuthal extent of the particles as pro- 
posed by Gold and Roelof (1976 ) for Jovian electrons. Quantitative 
studies are needed to determine if such an acceleration mechanism 
can produce enough particles to account for the increases seen at 
1 AU in 1973. A qualitative study of the energy requirements was 
discussed in Chapter VI. The energy flux in the solar wind was found 
to be ~ 2 X 10 6 larger than the energy flux in particles at 1.2 to 2. 4 
MeV streaming toward the sun, indicating the solar wind to be a pos- 
sible energy reservoir for interplanetary acceleration. 

3. Out of the Ecliptic. Another mechanism to create a posi- 
tive radial gradient is solar injection of protons at high solar latitude 
propagation out of the ecliptic to beyond 1 AU, and eventual diffusion 
back toward the sun in the ecliptic plane. Since little is known about 
out of the ecliptic processes, this possibility cannot be eliminated. 

4. Solar Flares. Another possible source for the increases 
seen by the Caltech experiment is solar flares, such as those which 
produce the prompt events but which, because of their different solar 
longitude, are not seen as prompt events. As discussed in Chapter 
X, the anisotropy observed late in prompt events is consistent with a 


diffusive streaming back toward the sun. Thus, if early in the prompt 
event the observer were poorly connected to the source of particles, 
the characteristic sharp rise with streaming from the sun might not 
be observed. Then if later the observer became better connected to 
the event due to rotation of the sun, for example, the event could be 
in its decay phase with the diffusive streaming back toward the sun. 

A simple model to produce this scenario divides particle- 
producing flares into two groups: those near 60° West solar longitude 
which are directly connected to the observer and produce the prompt 
events, and those to the east which do not become weii -connected 
until later in the event and produce delayed events. The model must 
be able to account for the relative number of prompt and delayed 
events. Twenty-five prompt events have been identified for the pe- 
riod included in the present study; they are listed in Table V-l. The 
diffusive anisotropy could be determined during the onset of 13 of 
these events. £ j| was negative, indicating flow from the sun, for 10 
of the 13 with values from -10 per cent to -100 per cent. For one 
period was consistent with 0 ; for the two remaining periods there 
was a large flow toward the sun (40 per cent and 100 per cent). 

Thus, for the majority of these events, the diffusive anisotropy is 
consistent with the solar source hypothesized. Eleven of the 25 
prompt events could be associated with a specific site of activity on 
the sun (Hurford, 1974). The site of the activity ranged from E57° 
to W80° solar longitude. Eight of the 11 were west of the central 

meridian. A typical maximum PEO rate for the 25 prompt events is 

2 

~ 0. 5 /second (i. e. , ~ 2/cm -sec-sr-MeV). 


The identification of delayed events is more subjective due to 
the lack of a sharp onset, so the determination of the number of such 
events is necessarily somewhat arbitrary. Let the threshold of a 
delayed event be a PLO rate of 0.01/second. The increases typically 
last several days, so that an estimate of the number of delayed events 
is the number of orbits from which at least one 6 -hour period is in- 
cluded in the final data set. Such periods are identified in Figure V-l. 
There are 26 orbits with at least one period. Another method is to 
count the number of transitions from PLO rates less than 0.01/sec- 
ond to PLO rates greater than 0.01/second which are not associated 
with identified prompt solar events. There are 36 such transitions 
from 72/273 - 74/2, The typical maximum PLO rate for these de- 
layed events is ~ 0. 05/second. Thus, there are somewhat more de- 
layed events than prompt events, and the prompt events are typically 
10 times larger than the delayed events. 

Whether solar flares can produce the delayed events depends 
on how much separation there must be from the flare site to the di- 
rectly connected longitude of ~ 60°W in order for the flare not to be 
seen as a prompt event. If the separation is comparable to the 
HWHM of 60° for the distribution in solar longitude of prompt solar 
events summarized by McCracken and Rao (1970 ) , then, to produce 
a comparable number of delayed and prompt events, the average po- 
sition of a delayed event will be ~ 60°E, It will take ~ 5 days before 
such an event becomes well-connected to earth, and by this time the 
event will have decayed to & 0. 01 times the size of a typical prompt 
event compared to the observed ratio of 0. 1 . 
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If, on the other hand, the required separation distance in order 
for an event not to be seen as a prompt event is only the 20° HWHM of 
the azimuthal spread of an event near 1 AU as deduced by McCracken, 
et al. ( 1971) , delayed events will be directly connected after ~ 2 
days and thus be nearly as large as prompt events. The current un- 
derstanding of the azimuthal propagation of prompt solar events does 
not indicate what the required separation is, so this solar flare model 
of delayed events remains as a possible explanation for the increases 
reported in this study. 

5. A Propagation Model . If the intensity increases observed 
in 1973 are due to particles being injected into the ecliptic plane at 
~3 AU either from out of the ecliptic or by local acceleration, their 
propagation back to I AU is closely analogous to galactic modulation. 

A numerical solution to the Fokker -Planck propagation equation in- 
cluding diffusion, convection, and adiabatic energy loss has been de- 
veloped using certain simplifying assumptions which are thought to be 
approximately correct. 

i) A steady state exists. Such a situation is suggested 
by the lack of dependence of £ u on shown in Figure 

VH-1. 

ii) Any injection is taking place beyond a distance L ~ 

3 AU such as in shocks or connection to an out of the ecliptic 
source. 

iii) Only the radial dependence is examined, and the ra- 
dial diffusion coefficient r is assumed independent of energy 

r r 

and radius. This dependence of x has been successfully used 
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Figure VII- 1 

The dependence of $, on (dU/dt)/U. (dU/dt)/U is approximated 
using finite differences in the PLO rate. Only 6-hour periods for 
which the PLO rate is known for the previous and subsequent 6-hour 
periods are included. The dashed line is the least squares fit to the 
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by Lupton and Stone (1973) to fit the temporal development 
of prompt solar particle events. 


Details of the numerical solution are presented in Appendix C. 
The solution has the form 

U(r) = Ae f|V/ *’ C ' r )Vr/* (?1) 


where f is a monotonic function of r varying from C at r = 0 to 1 at 

large r , Thus, the solution grows approximately exponentially with r 

with a scale length of ~ h/V. Using the typical solar wind velocity of 

440 km/second, C = 2. 77 (corresponding to y = -3. 15), and a H of 
2 1 2 

10 cm /second, the model produces a radial diffusive anisotropy 
and a modulation given by 

^DIF, x = 17 per cent, U(3 AU)/U(1 AU)= 11 . 

This diffusive anisotropy is comparable to the value of ~ 14 per cent 
found in this work. The modulation of 11 is consistent with the obser- 
vations of McDonald, et al. (1975 ). Smaller values of H can produce 
much larger modulation since the scale length is approximately pro- 
portional to h. These modulation factors can be reduced substantially 
by having the injection of particles occurring between 2 and 4 AU,and 
thus the gradient implied by eq. (7-1) would only hold out to 2 AU. An 
increasing x with radius would also reduce the modulation. 

The model is also consistent with other characteristics of the 
diffusive anisotropy found in this study. As shown in Appendix C, the 
diffusive anisotropy is approximately independent of x but proportional 
to V. H is not directly observable, but is thought to depend on the 
fluctuations in the magnetic field. No significant dependence of the dif- 
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fusive anisotropy on the fluctuations in the magnetic field was found. 

As shown in Figure VII -2, the diffusive anisotropy did increase with 
increasing solar wind velocity, approximately as predicted by the 
model. 

Thus , there is a simple model for propagating particles in- 
jected at ~ 3 AU which explains the size of the diffusive anisotropy 
found at 1 AU, produces a radial gradient such as those indicated by 
the measurements of McDonald, et al. (1975 ) for a reasonable size of 
the diffusion coefficient, and whose predictions of the dependence of 
the diffusive anisotropy on H and the solar wind speed are consistent 
with the observations. 

6. Role of Continuous Solar Injection. The new measurements 
reported herein indicate that continuous solar injection is not neces- 
sary to produce increases in the low-energy proton flux such as those 
seen during 1973. The strongest direct evidence that some events are 
due to continuous solar injection is found in the observations of Krim- 
igis, et al. (1971) and the detailed analysis of these data by Roelof 
and Krimigis (1973). The observations covered 3 solar rotations in 
1967. Three delayed events are discussed in which there is a large 
anisotropy throughout the event due to particles coming from the sun, 
indicating continuous solar inj ection. These data are at energies 
£ 0. 3 MeV. 

At energies nearer the 1.3 -2.3 MeV used in the present 
study, the direct evidence consists of anisotropy measurements made 
by Fan, et al. (1968) using two periods, each of ~ 6 hours. The 0.8 - 
2. 0 MeV protons were found to be streaming from the sun in a direc- 


Figure VI 1-2 

The average x-component of the diffusive anisotropy as a function of 
the solar wind speed. The solid curves are calculated using the prop- 
agation model discussed in Section VII-A-5 for the parameters indicated. 
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tion close to the long-term average magnetic field direction. These 
data indicate occasional solar injection, but do not require continuous 
injection. 

Consequently, at least at energies £ 1 MeV, continuous solar 
injection may not be the predominant source over the solar cycle of 
the proton increases seen between prompt solar events; rather, some 
other mechanism, such as those discussed above to explain the in- 
creases seen in 1973, may be the usual source. 

B« Quiet-Time Measurements of Rao, et al, (1967b) 

The anisotropy measurements of Rao, et al. (1967b) were taken 
during the periods between observed increases. About 66 days of data 
from the UTD experiments on Pioneer 6 and 7 during 1965 and 1966 
were used. As noted in Chapter I, the flux during these periods, when 
extrapolated to the energies used in the current study, is within the 
range of fluxes used in the present study. In contrast to the typical 
observed anisotropy of 7 per cent reported in Chapter VI, Rao, et al. 
reported a;i anisotropy of 0. 19 $ ± 0. 05 $ . Forman and Glees on 

(1975 ) have discussed the difficulties in reconciling this small aniso- 
tropy with current understanding of propagation. A possible explana- 
tion of this result is that the observations included a large isotropic 
background. This possibility is consistent with the fact that the 
smallest flux seen by the University of Chicago experiment (Fan, et 
al. , 1968 ) at the same energies during 1966 was ~ 30 times smaller 

than that seen by Rao, et al. Even lower fluxes have been seen at 
these energies by the Caltech EIS experiment on IMP -7 during 1973 
(Mewaldt, et al. , 1975b). This background in the UTD instruments 
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will minimally affect anisotropy measurements reported by the UTD 
group at higher fluxes such as those seen during prompt solar particle 
events. 

C. Magnetic Field Direction 

The interplanetary magnetic field has long been thought to play 
a dominant role in low-energy proton propagation. McCracken, et al. 

( 1968 ) investigated the dependence of the observed anisotropy on the 
magnetic field direction using 1-hour averages during the early part 
of prompt solar particle events. A strong dependence was found, and 
was interpreted by the authors as consistent with an observed aniso- 
tropy consisting of two parts: an invariant convective part and a field- 
aligned diffusive part. This picture of the diffusive anisotropy being 
field-aligned has been used by many authors, e. g. McCracken, et al. 
(1971), McKibben (1973 ), Wibberenz (1974). Rec ?nt observations 
have raised questions about this simple picture. Allum, et al. (1974) 
found that during the easterly anisotropy observed late in the decay of 
solar particle events, the observed anisotropy direction was independ- 
ent of the observed magnetic field direction. Six-hour averages were 
used. The authors noted that this unexpected result could be ex- 
plained in terms of field -aligned diffusion if the diffusion coefficient 

along the magnetic field direction had the dependence: 

2 _ 

K || « 1/cos (4>-4> ) 

where 4> is the long-term average of the magnetic field direction. 

No justification for this dependence was presented. Pesses and 
S arris (1975 ) have reported periods when they deduced the diffusive 
anisotropy could not be field -aligned, using 15-minute averages by 
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examining periods when the magnetic field was radial. 

The previous observations have used the direction of the ob- 
served anisotropy which includes the convective anisotropy, which is 
unrelated to the magnetic field direction. For a direct determination 
of the dependence of the diffusive anisotropy direction on the magnetic 
field direction, | CON mast be subtracted from £ DIF . as was done in 
the present study. 

If the diffusive anisotropy were field -aligned, the following 
functional dependence is expected: 

h = a + b 4> , 

*DIF 

with a consistent with 0 and b consistent with 1 . If diffusion were 
isotropic, b would be consistent with 0. 

The best least-squares fit to the data was found for 
a = -10. 2° ± 2. 9° , b = 0.71±0.06 . 

Thus, a strong dependence of the diffusive anisotropy direc- 
tion on the magnetic field direction has been found, indicating K x to 
be typically less than Hu for these 212 periods. The difference from 
the dependence expected for field-aligned diffusion is consistent with 
a finite kJk\\ . The size of this ratio depends on the relative sizes of 
the spatial gradients perpendicular and parallel to the magnetic field. 
Since information from only one satellite is available, these gradients 


are not known. 
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APPENDIX A. ANISOTROPY MEASUREMENTS 

1, Introduction 

This appendix describes the determination of anisotropies 
used in this study and discusses the statistical uncertainties of such 
measurements. The IMP-7 instrument accumulates counts separately 
for each of 8 sectors. An anisotropy is calculated using the accumu- 
lated counts, and the statistical uncertainties associated with the ac- 
cumulated counts produce statistical uncertainties in the computed 
anisotropy. Ideal instruments are assumed which require no cor- 
rections for dead times. Roelof (1974 ) has discussed dead time 
effects. The correction for finite opening angles will be discussed in 
Section A-4. Only first order anisotropies are considered. 

2. Calculation of the Observed Anisotropy 

A cosine expansion is fit to the accumulated counts for the 

different sectors. Let Y. be the number of counts for sector i . 0. 

l i 

the statistical uncertainty associated with Y^ , and <}>♦ he the average 
angle of sector i The angles 4>.. for EMP-7 are given in Figure II-2. 
Then the Y^'s can be approximated by the function 

f x (fy A,£> 4>) = A(l + 6cos(4>.-4>)) . (A- 1) 

Following the method of Zwickl and Webber (1974), a least-squares 
fit is made to an equivalent function 

f-(4>., A, tt, (3) = A + CXcos <j>. + (3 sin 4>. . (A-2) 

C* \ 11 

x 

The least-squares fit has been done using a . = (Y..) 2 . The parame- 
ters of eqs. (A-l) and (A-2) are related by 

i = (a 2 +p 2 )*/A 

4? = tan"*(p/tt) . 


(A-3) 

(A-4) 
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Although weighted averages have been used to determine A, a, and 
(3 , the approximation that all CT's are equal has been used to deter- 
mine the uncertainty in the fit parameters. The approximation is a 
good one for small anisotropies. Using this approximation, the best 
fit is made for 


A = (l/s)£ Y. , where s is the number of sectors 

a = (2/s )£ Y i cos 

P = (2/s)S Y.sin<j>. 

The observed anisotropy amplitude is 

f _ = w£ , where w = w w 

9 OBS b s a 

The smoothing factor w , which equals (ir/s )/sin(Tr/s) , corrects for 

s 

the reduction in the anisotropy due to the finite number of sectors 
(Chapman and Bartels, 1940 ). w a is the smoothing factor correct- 
ing for the finite opening angle of the instrument, and is discussed in 
Section A-4. 

3. Statistical Significance 

To understand the significance of an observed anisotropy, the 
likelihood of possible true anisotropies must be determined. The 
probability of the true anisotropy, given an observed anisotropy, is 


determined using Bayes* Theorem: 

P <W'*obs > = ,lP (*trae ),p( *obs'*true ) (A ' 5) 
where v is normalization. It is assumed that the true anisotropy 

distribution is of the form given by eq. (A-l)or equivalently eq. (A-2). 

Certain approximations are made which are valid only for small 

anisotropy amplitudes. 

Let the true distribution be given by eq. (A-2) and character- 


135 


ized by A , a , and 3 . Let 

* n n n 


i 


i 

2 2 z 

<■»+?.) 


true 


o * o 

~x — 


and rf> = tan"^(3 /a ). The distributions of Y.'s sampled from this 
y o r o o i r 

distribution are approximated as Gaussians, each having an equal 
standard deviation 0 . For the present work, <7 = (s/E l/o^ ) 2 . 
Each sample distribution is characterized by A, a, and (3 . The ap- 
proximation is made that the relative uncertainty in A is negligibly 
small compared to the uncertainties in a and (3 , which will be the 
case for | « 1. a and (3 then have Gaussian distributions with 


means of a Q and (3 q respectively, and standard deviations of a( 2/s) 

Define a scale factor for the anisotropy amplitude: 

i 

z = w (2/s) 2 cr/A Q . 

Let 


r = W z * 

x = hrae /z ' 

z is approximately the anisotropy amplitude expected to be measured 
if the true distribution is isotropic. Hence, r is a measure of the sta- 
tistical significance of an observation -- the larger r , the more sig- 
nificant the observation. 

Equation (A- 5) becomes 

P(x, $ o |r, 4>) = k^fx, <J> Q )P(r, $|x, <J> q ) , (A-6) 

and from the Gaussian distributions at a and p , 

2 2 

P(r, <J>|x, ^> o ) = (r/2ir) e *V" . (A-7) 

Integrating over angle gives the probability of measuring r given a 
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true distribution characterized by x: 

i 2 l 2 , 

F(r|x) = re -2 ** e“ z3C I Q (xr) , (A-8) 

where I is the zero -order modified Bessel function. The mode and 
o 

68. 2 per cent and 95. 4 per cent confidence intervals of this proba- 
bility function are plotted in Figure A-l as a function of x. There is 
a 68. 2 per cent (95.4 per cent) probability that an observed r will be 
within the 68.2 per cent (95.4 per cent) confidence interval. The in- 
tervals are chosen to minimize their lengths. Note that as x -* 0 
(i. e. , small true anisotropies) the most likely measured anisotropy 
approaches z , not 0 . 

In order to compute P(x, (J> Q |r, , an assumption must be 

made about the a priori probability of different anisotropies, P(x, 4> q ). 
The assumption has been made that P(x, 4> q ) is independent of x and 
. Other assumptions are possible, however. Another assumption 
might be that all points in the a - P plane are equally likely, i. e, , 

P(x, <J> Q ) ~ x. This a priori assumption would be appropriate if the 
two orthogonal directions corresponded to two independent sources of 
particles. The assumption which is used of a priori equally likely 
amplitudes corresponds to a single particle source in an arbitrary di- 
rection, and so is considered more appropriate to inter pi erne tar y 
propagation. 

So now 

P(x, 4> 0 |r, cj>) = k 2 P(r, <|>|x, <f> Q ) . (A-9) 

Integrating eq* (A-9) over <f> Q gives the probability that the 
true distribution has an anisotropy amplitude of x given an observa- 


tion of r : 


137 


Figure A-l 

The 68.2% and 95.4% confidence Intervals for observing an anisotropy 
amplitude given a true anisotropy amplitude. Both the observed and 
true anisotropy amplitudes have been normalized by z, where 

z = w (2/s) 0 ’ 5 o/A 0 . 

w is the smoothing factor, usually =1., which depends on the number of 
sectors s and the angular response of the detector system, a Is the 
statistical uncertainty in the number of counts in a sector (assumed 
to be Independent of sector). A Q Is the average number of counts per 


sector. 
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_i 2 . . ^ 2 _ 

P(x|r) = (2/ir) 2 e~ r e~ aX I o (xr)/I Q (r /4) . (A-10) 

This probability distribution is plotted for r = 0, 2, and 3 in Figure 
A -2. The 68. 2 per cent and 95. 4 per cent confidence intervals are 
also indicated. The most likely x and the 68. 2 per cent and 95. 4 per 
cent confidence intervals are plotted as a function of r in Figure A -3. 
For r £ tJ'T the most likely true anisotropy amplitude is 0. For 
r £ 3 the most likely true amplitude is nearly the observed amplitude 
with a standard deviation of nearly z . 

Integrating eq. (A -9) over the anisotropy amplitude x gives 
probability of the true direction being d D given an observation of r 
and d> : 

2 sin^ Ad 

P(«> 0 |r.4>) = k 3 e [l + erf(- r - ^ 8 -A* )] (A-ll) 

where Ad = d -d * and k is a normalization constant. The distribu- 
o 3 

tion is symmetric with respect to Ad and the most likely Ad is 0. 

For large r , that is, for statistically significant observations, 

i 2 2 

P(4> o |r, «j>) a k 3 e‘ 2r A<l> (A-12) 

so that the probability distribution for the true angle approaches a 
Gaussian distribution with a standard deviation in radians of 1/r = 

z ^OBS * 

Figure A -4 shows the half -width - half-maximum for eq. (A- 

11) as a function of 1/r . For r < 0.431 the HWHM is not defined 

because the distribution is too broad. The figure shows that the 
1 

HWHM of (2 07i2) 2 /r appropriate for large r [eq. (A-12)] is a good 
approximation for r as small as 0.5. 


Figure A- 2 

The probability distribution for the normalized true anisotropy ampli- 
tude for three values of the normalized observed anisotropy amplitude. 
Also shown are the 68.2% and 95.4% confidence Intervals for each dlstrl 
butlon. The normalization factor z Is defined In the caption of Fig- 


ure A-l. 
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Figure A-3 

The 68,2% and 95.4% confidence Intervals for the normalized true anisot- 
ropy amplitude given a normalized observed anisotropy amplitude. The 
normalization factor z Is defined In the caption of Figure A-l. 
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Figure A-4 

The half-wldth-half-maxlmum of the probability distribution of the 
difference In the direction of the observed anisotropy and the direct- 
ion of the true anisotropy as a function of the normalized observed 
anisotropy amplitude, r. The HWHM cannot be defined for r < 0.43. 
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4. Finite Opening Angle 

Sentman and Baker (1974) have shown that the observed 
counting rate is given by 

C(v) = 4irS 2 -jj- D^S^tcosv) , 

where are Legendre polynomials, are the coefficients in Le- 
gendre polynomial expansion of the pitch angle distribution, S . are 

As 

the coefficients in the expansion of the detector angular response, and 
V the angle between the direction of the anisotropy and the symme- 
try axis of the detector teles "ope. 

The first order anisotropy (l = 1) is reduced by 

7T 

T S(0) cos 8 sin 6 d8 

— - — • cos(0 ) , 

w Tf 'o' * 

J S{0) sin0 d0 

where 0 q is the elevation of the direction of the anisotropy above the 
plane of rotation cf the telescope. For the Caltech instrument this 
plane is.tfte ecliptic plane. 0 q is assumed to be 0 for this study, so 
the measured anisotropies are the projection of the true anisotropy 
vector onto the ecliptic plane. 

S(8) has been determined for the Caltech EIS experiment on 
IMP-7 using Monte-Carlo simulation, w is found to be 1.0/0. 951 . 
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APPENDIX B. DETERMINATION OF ENERGY SPECTRA 
1. Method 

The procedure used for calculating proton energy spectra has 
been discussed by Hurford (1974). This appendix discusses the par- 
ticulars of the procedure used in this work. 

The range of D2 pulse heights for protons stopping in D2 has 
been divided into 9 bins and the number of analyzed events with the 
corresponding pulse heights accumulated for each of the 6 -hour peri- 
ods used in this study. Corrections, which are discussed below, are 
made to the accumulated counts: 

N. = C.- A. -f.(P25.+P256.) (B-l) 

where N. is the corrected number of counts for the i th bin C is 
i * i 

the observed number of counts, A. is the alpha contamination, f. is 

1 i 

the foldback correction factor, and P25 i and P256 i are the number 
of particles which penetrate D2 and trigger D5 and D56 respectively, 
and have a D2 pulse height corresponding to the i th bin. The energy 
interval of incident protons for each of the pulse height bins is deter- 
mined, and a least-squares fit is made to the function 

dN/dT = AT y , 

where T is the incident kinetic energy of the proton. 

During the period 73/86 to 73/154 the offset of the D2 pulse- 
height analyzer differed from the nominal value by — 4 channels -- 
the nominal channel 16 becoming channel 20. This offset was deter- 
mined using the D2 pulse height for D5 single events. The nominal 
electronic calibrations given by Mewaldt and Vidor (1976) have been 
adjusted by 4 channels for this period. In addition, the lowest two D2 


pulse height bins have not been used in order to maintain the same 
incident energy threshold. 

2. Incident Energy Intervals 

The D2 pulse heights and the corresponding incident particle 
energies are listed for the 9 proton bins and 3 alpha bins {see below) 
in Table B-l. The incident energy is determined by adding the energy 
loss in D2 to the calculated energy loss in the thin aluminized mylar 
window covering the detector stack. The energy loss in D2 is deter- 
mined from the D2 pulse height using the results of the electronic 
calibration of the instrument {Mewaldt and Vidor, 1976). The 
equivalent thickness of the aluminized mylar window is 2.4 mg/cm 
of lucite as measured using 6.051 MeV and 8.785 MeV a-particles 
from 2 l 2 p b * Since (sec 0> , where 0 is the angle of incidence with 
respect to the symmetry axis of the detector stack, is 1.054 as de- 
termined by a Monte -Carlo simulation of the detector stack response, 

the effective thickness of the window for particles stopping in D2 is 
2 

2. 53 mg/cm . From this effective thickness and the known residual 
energy loss in D2, the incident energy is calculated using the range - 
energy tables by Janni (1966). Alpha energies are determined using 
the relationship: 

RJT) = (m /4)R (T/m > (B-2) 

a a p a 

where R a is the alpha range, R^ the proton range, T the alpha 

kinetic energy, and m & the alpha mass. 

Errors in the determination of the effective window thickness 

produce errors in the determination of v and hence in Sta- 

3 CON 

tistical errors in the measurement of the thickness or determination 

















i 
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TABLE B- 1 


Bin No. 

D2 Pulse 
Height 


Incident Energies 
MeV /nucleon 

protons 

channels 

nominal 

73/86 - 73/154 

1 

16 - 16 

1. 185 - 1.204 

not used 

2 

17 - 19 

1.204 - 1. 261 

not used 

3 

20 - 24 

1. 261 - 1. 379 

1. 185 - 1.282 

4 

25 - 29 

1. 379 - 1. 515 

1.282 . 1.406 

5 

30 - 34 

1.515 - 1.667 

1. 406 - 1. 544 

6 

35 - 39 

1. 667 - 1. 829 

1. 544 - 1. 698 

7 

40 - 44 

1. 829 - 1. 999 

1. 698 - 1. 861 

8 

45 - 49 

1. 999 - 2. 172 

1. 861 - 2. 031 

9 

alphas 

50 - 53 

2. 172 , 2. 315 

2. 031 - 2. 172 

10 

87-119 

1.44 - 1. 69 

1.41 - 1.66 

11 

120 - 149 

1. 69 - 1. 93 

1.66 - 1.91 

12 

150 - 198 

1. 93 - 2. 37 

1. 91 - 2. 34 



TABLE B-2. 

Alpha Energies Which Contaminate Proton Bins 

Bin No. 

Nominal 

73/86 

- 73/154 


Mean 

Width 

Mean 

Width 


, MeV 

. , MeV . 

/ MeV 

( MeV . 


nucxeon nucleon 

' nucleon ’ 

* nucleon' 

1 

1. 081 

0. 002 

not used 


2 

1.086 

0. 008 

not used 


3 

1.097 

0.014 

1.086 

0. 012 

4 

1. 113 

0.019 

1. 100 

0. 016 

5 

1. 132 

0.019 

1. 116 

0.016 

6 

1. 153 

0.022 

1. 132 

0.020 

7 

1. 176 

0.024 

1. 158 

0.023 

8 

1.201 

0. 026 

1. 181 

0.024 

9 

1. 223 

0.019 

1.204 

0.021 
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of (sec 6) are quite small (~0.2 per cent). However, there may be 

systematic errors in the range -energy table, in the approximation of 

mylar as lucite, or in eq. (B-2). An estimate of these systematic 

errors is made by using the range -energy tables of Northcliffe and 

Schilling (1970) for alpha particles in mylar to determine the window 

2 

thickness. A thickness of 2.47 mg/cm is determined, compared to 
2 

2.4 mg/cm for the previous calculation. If the thickness were 2.47 
instead of the 2. 4 used, a \ of -3.00 would be measured as a y of 
-3.04, producing an error in on ^ ~ 0* 2 per cent. Conse- 

quently, possible errors in the determination of the window thick- 
ness have at most negligible effect on the anisotropies of this study. 

3. Alpha Correction 

Some of the counts in the 9 proton bins are due to stopping 
alphas. An estimate of this contamination is provided by the number 
of particles that lose too much energy in D2 to be stopping protons. 

A least-squares fit of the form 
dN/dT = AT V 

where T , alpha incident energy, is in MeV/nucleon, is made to the 
three alpha bins shown in Table B-l. This alpha spectrum is then 
extrapolated to lower energies to compute the number of counts that 
are subtracted from the observed counts in the 9 proton bins, which 
are designated A.. in eq. (B-l). The average alpha incident energy 
and bin width corresponding to each of the 9 proton bins are shown in 
Table B-2. 

4. Foldback Correction 


Some of the particles contributing to the 9 proton bins pene- 


r 
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TABLE B-3. D5 Pulse Heights Used to Bin D256 Events 


Bin No. 

D5 Pulse Height 

73/86 - 73/154 


Nominal 


1 

234 - 400 

not used 

2 

165 - 233 

not used 

3 

102 - 164 

149 - 400 

4 

65 - 101 

89 - 148 

5 

44 - 64 

61 - 88 

6 

31 - 43 

42 - 60 

7 

19 - 30 

28 - 41 

8 

10 - 18 

17 - 27 

9 

3 - 9 

7 - 16 


TABLE B-4. Foldback Correction Factors (L ) 


Bin No. 

Nominal 

73/86 - 73/154 

1 

0. 10 

not used 

2 

0.20 

not used 

3 

0. 36 

0.20 

4 

0.87 

0.46 

5 

1. 50 

0.98 

6 

2. 05 

1. 60 

7 

2.20 

2.08 

8 

1. 81 

2. 13 

9 

2.90 

1. 90 
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trate D2 and stop in the thin dead layers on the inside of the annular s 
D3 and D4. The size of this correction depends on the spectrum of 
penetrating particles and the efficiency with which penetrating parti- 
cles are rejected by D3, D4, or D5. 

The spectrum of penetrating particles is determined by accu- 
mulating P25^ -- the number of D25 events whose D2 pulse height 
corresponds to one of the 9 proton bins listed in Table B-l. Due to 
an instrumental anomaly, D25 events are occasionally read out as 
D256 events with the D2 pulse height lost. Such events are recognized 
by their zero D6 pulse height. Such events have been binned according 
to the average D2 pulse height appropriate for the observed D5 pulse 
height, as shown in Table B-3. The accumulated counts are desig- 
nated P256. . 

1 

The sum of P25^ and P256. is then multiplied by a foldback cor- 
rection factor f^ which is the ratio of the geometric factor of D2 5 
events to the effective geometric factor of particles stopping in the 
dead layer of D3 or D4. The efficiency of D3 or D4 in rejecting par- 
ticles depends on the particle's energy, so there is a different f. for 
each of the 9 bins. The f^ were determined by comparing the energy 
spectrum of the Caltech experiments on IMP -7 and IMP -8 during 0 UT 
to 1200 UT 75/234. This period was chosen because both IMP-7 and 
IMP-8 were outside the magnetosphere and the relative number of 
penetrating to stopping particles in D2 was large. The energy spec- 
trum observed by IMP -8, which needs no foldback correction, was 
well fit by a power law. The size of the foldback correction for IMF- 
7 was chosen to make the corrected spectrum equal to the extrapola- 


153 


V IMP 7 D25 power law energy spectrum. The T 

tion of the observed IMP-7 u“P 

, n Table b . 4. The uncertainty of the T is estimated by 
are given m Table a 

• ^ fKe above period to the values 

comparing the values determined using the above p 

• a 12-00 to 24:00 UT 75/235. The two sets of f. agreed 

determined using 12.00 to 

wiihia ~ c„. e~.pt <9 ™ ~ M ' • 

,, and , 2 , which diitered bp a «- •< - h. «* h.,. c^bi. 

effect on the spectra. 

5. Example 

— nr0 cedure described in this ap- 

The results of apply in g tlie P 

„ f from 73/345 are shown in Figure B-l. The uncorrected 

pendix to data from 73/5«= a 

.p.c„». » “* ““ ‘ 

cent— ibib— -ion b— « >•«“>" “« ^ “ P, “ 

Ration i. comparable .. CM*-*-- 

rection is somewhat smaller than typical, as seen by the y o - . 

comoared to mean y of -3. 15. 




yrrV 




Figure B^l 

The determination of the spectral index y for 73/345. The observed D2 
counts, the correction due to o-contami nation, the correction due to 
the foldbaek of protons which penetrate D2, and the corrected D2 counts 
are shown as a function of energy. 
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APPENDIX C. PROPAGATION MODEL 


The propagation of cosmic rays in interplanetary space is de- 
scribed by the Fokker -Planck equation: 


= VM** VU)-V*(UV)+±V*v(^ (a<T)TU)) . (C-l) 

This equation includes the effects of diffusion, convection, and adia- 
batic energy change. The equation is solved for r e[0, L 1 using cer- 
tain simplifying assumptions: 

1) all variables and parameters are independent of time; 

2) the solar wind velocity V is radial and independent of position; 

3) the energy spectrum of U is a power law; 

4) a = 2 ; 

5) U remains finite as r -* 0 ; 

6) no source of particles for r < L ; 

7) ^ is diagonal in reference frame aligned with the radial direc- 
tion; Xq = x^ = 0 and n rr = h is independent of energy and 
position. 

Let 

U(r, 6, $,T) = R(r)Q(6, 4>)T V '* . (C-2) 

Equation (C-l) separates and R is given by 


R" + (- - ~)R ! - R = 0 . 
v r K r xr 


Let 


Then 


R(r) = y(r)e 


Vr/H . 


y" + y'(g± .YLV) . y ,(.?-! )2^V/K)r = Q . 


(C-3) 

(C-4) 

(C-5) 


This equation is solved using power series solution, 
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= r 8 S a r 1 
n=0 n 


(C-6) 


The indicial equation for s has roots s = 0, -1 , so one solution is ob- 


tained of the form 


uu 

yi (r) = S a r n . 
n=0 


(C-7) 


Substitution of eq. (C-7) into (C-5) determines the coefficients a 

1 1 


a j = a o (C-l)V/H 


a 2 = a l /< 3a 0 ) 


, • £o. n = 3 

n +n 


(C-8) 


where a Q is arbitrary. 

An independent solution exists with form 


i 00 

y 2 ( r ) = Ay^( r)Sn r + r" 1 Zj b r n 

n=0 r 


y 2 is infinite at the origin and so is discarded. Thus, the desired 


solution is given by 


R(r) = e Vr/x S 


a o is determined by the value of R(L) . 

Physically, the density at r = L is set by the source of parti- 
cles at r a L. The solution then shows the radial distribution that 


will be established by the propagation of these particles. 

The solution, normalized to 1 at r = 0 , is plotted for several 


values of the Compton -Getting factor in Figure C-l. The radial dis- 


tance has been scaled using x = r/(x/V). For x. = 10 21 cm 2 /sec and 


V = 440 km/sec, k/V = 1.5 AU. 


ir; ^ -I! r*V 



1 





) 
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Figure C-1 

The radial solution as a function of x for four values of the spectral 
index y. The solutions are normalized to a value of 1.0 at x * 0. 

The Compton ^Getting factor, C, corresponding to each y is given. 



RADIAL SOLUTION 

fC = 4.5 

/ / ly«-5.75 

f C = 2.767 
/ 1 X= -3.15 




» 3 h f dR 


’DIF 


w 


(#y 


R = 


160 

-3V dkR 


w 


"3x 


(C-9) 


so that the slope of a solution plotted in Figure C-l is proportional to 

^ , * -3CV 

the diffusive streaming. At x = 0 , ^ ■ - = C ana - — • 

d0nR - . -3V 


For large x , 


-* 1 and £. 


~d > c * *DIF w 

The dependence of | DJF computed using eq. (C-9) on V and x 

is shown in Figure C-2. The anisotropy is computed at 1 AU for 
1. 675 MeV protons with y = -3. 15 . These values are typical of the 
data used in this work. S DIF is relatively independent of x , in- 
creasing by less than a factor of 2 when K increases by a factor of 
40. In contrast, | is nearly proportional to V . A qualitative 
picture is that the diffusive streaming approximately offsets the con- 
vective streaming in this steady-state model. Hence, increasing V 
and thus the convective streaming also increases the diffusive stream- 
ing. On the other hand, increasing x does not increase the convective 
streaming so the radial gradient adjusts to keep the diffusive stream- 
ing approximately unchanged. 


«vPRO0TJCrolLITY0FTEa 


Figure C-2 

The dependence of calculated diffusive anisotropy amplitude 5 on the 

radial diffusion coefficient < rr and the solar wind speed V sw . For the 

upper panel V sw is set and K rr varied. The asymptotic limits on £ for 

large and small < are indicated. For the lower panel k is set and 
rr r rr 

Vsw vaned * 




APPENDIX D. DATA SET 

Table D-l lists the 6 -hoar periods ased for the results of this 
study. The criteria used for their selection have been discussed in 
Chapter V. Each period is identified by yy ddd q , where yy is the 
year, ddd is the day number, and q is a number from 0 to 3 which 
specifies one of the four 6 -hour periods of the day. In addition, the 
FLO rate in counts /second is given for each period. 
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TABLE D-l. Periods Used 


TIM8 PIC 


72 

285 

0 

0.0088 

73 

72 

285 

l 

0.0150 

73 

72 

285 

2 

0.0154 

73 

72 

285 

3 

0.0153 

73 

72 

294 

0 

0.0045 

73 

72 

294 

1 

0.0048 

73 

72 

294 

3 

0. 0036 

73 

72 

317 

2 

0.0232 

73 

72 

317 

3 

0.0205 

73 

72 

322 

2 

0.0031 

73 

72 

342 

2 

0.0044 

73 

72 

344 

0 

0.0031 

73 

72 

345 

1 

0.0039 

73 

72 

345 

3 

0.0051 

73 

72 

346 

0 

0.0067 

73 

72 

346 

I 

0.0151 

73 

72 

346 

2 

0.0059 

73 

72 

346 

3 

0,0054 

73 

72 

355 

1 

0.0068 

73 

72 

355 

2 

0.0058 

73 

72 

355 

3 

0.0043 

73 

72 

356 

0 

0.0037 

73 

72 

356 

l 

0. 0040 

73 

72 

356 

2 

0.0049 

73 

72 

356 

3 

0.0051 

73 

72 

357 

0 

0,0041 

73 

72 

357 

1 

0.0057 

73 

72 

357 

2 

0.0086 

73 

72 

357 

3 

0.0054 

73 

72 

359 

2 

0.0040 

73 

72 

359 

3 

0.0100 

73 

72 

360 

0 

0.0085 

73 

72 

360 

1 

0.0082 

73 

72 

360 

2 

C, 0298 

73 

72 

360 

3 

0.0234 

73 

72 

361 

0 

0.0174 

73 

72 

361 

1 

0.0140 

73 

72 

361 

2 

0.0143 

73 

73 

2 

2 

0.0031 

73 

73 

2 

3 

0.0052 

73 





PLO 


TIMP 


3 

0 

0, 

.0044 

73 

53 

3 

3 

1 

0, 

.0057 

73 

54 

0 

3 

2 

0< 

.0072 

73 

92 

0 

3 

3 

0, 

.0103 

73 

92 

1 

4 

0 

0, 

.0163 

73 

92 

2 

4 

l 

0, 

.0206 

73 

93 

0 

4 

2 

0, 

.0313 

73 

93 

1 

4 

3 

0, 

. 0332 

73 

93 

2 

5 

0 

0. 

.0265 

73 

93 

3 

5 

1 

0, 

,0189 

73 

94 

0 

5 

2 

0, 

► 0139 

73 

94 

1 

5 

3 

0, 

,0137 

73 

94 

3 

6 

0 

0. 

► 0115 

73 

95 

0 

6 

1 

0. 

► 0121 

73 

95 

l 

6 

2 

0. 

0224 

73 

117 

0 

6 

3 

0. 

, 0339 

73 

117 

2 

7 

0 

0. 

, 0547 

73 

117 

3 

7 

1 

0. 

C892 

73 

128 

1 

T 

1 

2 

0. 

0968 

73 

128 

2 

7 

3 

0. 

0623 

73 

128 

3 

8 

0 

0. 

0678 

73 

129 

0 

19 

3 

0. 

0051 

73 

129 

1 

20 

0 

0. 

0056 

73 

129 

3 

20 

1 

0. 

0102 

73 

130 

0 

20 

2 

0. 

0052 

73 

130 

l 

20 

3 

0. 

0070 

73 

130 

2 

31 

0 

0. 

0036 

73 

130 

3 

31 

1 

0. 

0072 

73 

131 

0 

31 

2 

0. 

0171 

73 

131 

l 

31 

3 

0. 

Cl 97 

73 

131 

*> 

32 

0 

0. 

0274 

73 

180. 

0 

32 

1 

0. 

0172 

73 

i32 

1 

32 

2 

Cc* 

C093 

73 

132 

2 

32 

3 

0. 

C190 

73 

133 

0 

33 

0 

0. 

0155 

73 

133 

1 

33 

1 

9. 

0307 

73 

133 

2 

33 

2 

0 . 

0344 

73 

141 

0 

53 

0 

0 . 

0686 

73 

141 

2 

53 

1 

Oo 

0518 

73 

141 

3 

53 

2 

0 . 

0236 

73 

143 

3 


PLO 

0.0X63 
0.0157 
0.0243 
C.02C6 
0.0166 
0.0096 
0.0121 
0.0122 
0.0089 
0.0057 
0.0039 
0.0061 
0.0052 
0.0072 
0.0079 
0.0033 
0.0042 
0.1102 
0.0951 
0. 0781 
0.0888 
0.0680 
0.2141 
0.1677 
0.1089 
0.0747 
0.0535 
0.0611 
0 . 06 86 
C. 0623 
0.0742 
0.0920 
0.1350 
0.2553 
0.2424 
0. 16C4 
0.1482 
C. 0261 
0.0165 
0.0277 


r 






1 



J 

t * 


t 
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TABLE D- 1. Periods Used (continued) 



TIME 


PLO 


TIME 


73 

144 

0 

0.0312 

73 

221 

0 

73 

144 

1 

0.0238 

73 

221 

1 

73 

144 

3 

0.0235 

73 

221 

2 

73 

145 

0 

0.0062 

73 

221 

3 

73 

145 

1 

0.0081 

73 

222 

0 

73 

145 

2 

0.0076 

73 

222 

2 

73 

145 

3 

0.0071 

73 

222 

3 

73 

146 

0 

0.0085 

73 

223 

0 

73 

146 

I 

0.0076 

73 

229 

t 

73 

146 

2 

0.0121 

73 

229 

2 

73 

146 

3 

0.0204 

73 

229 

3 

73 

157 

0 

0.0427 

73 

230 

0 

73 

159 

1 

0.0254 

73 

230 

2 

73 

159 

2 

0.0224 

73 

230 

3 

73 

166 

3 

0.2044 

73 

231 

0 

73 

167 

0 

0.0952 

73 

231 

2 

73 

167 

1 

0.0690 

73 

233 

2 

73 

167 

3 

0.0612 

73 

234 

3 

73 

168 

0 

0.0476 

73 

235 

0 

73 

168 

1 

0.038! 

73 

235 

1 

73 

168 

2 

0.0486 

73 

235 

2 

73 

168 

3 

0.0388 

73 

242 

1 

73 

169 

0 

0.0214 

73 

242 

2 

73 

169 

l 

0.0050 

73 

242 

3 

73 

192 

0 

0.0210 

73 

243 

0 

73 

192 

1 

0.0386 

73 

243 

1 

73 

192 

2 

0.0553 

73 

243 

2 

73 

192 

3 

0.0579 

73 

243 

3 

73 

193 

0 

0.0710 

73 

244 

0 

73 

193 

1 

0.0997 

73 

244 

l 

73 

193 

2 

0.1209 

73 

244 

2 

73 

193 

3 

0. 1406 

73 

244 

3 

73 

194 

0 

0.1538 

73 

245 

0 

73 

196 

1 

0.0560 

73 

245 

1 

73 

219 

2 

0.0331 

73 

245 

2 

73 

219 

3 

0.0269 

73 

245 

3 

73 

220 

0 

0.0237 

73 

247 

0 

73 

220 

1 

0.0243 

73 

247 

2 

73 

220 

2 

0.0366 

73 

247 

3 

73 

220 

3 

0.0350 

73 

248 

0 


PLO 


TIME 


PLO 

0.C291 

73 

240 

1 

0.0281 

0.0266 

73 

256 

1 

0.2961 

0.0184 

73 

256 

2 

0.2774 

0.0186 

73 

257 

2 

0.2828 

0.0202 

73 

2 57 

3 

0.2433 

0.0179 

73 

258 

0 

0.1793 

0.0151 

73 

258 

1 

0.1417 

0.0162 

73 

258 

2 

0.0922 

0.0175 

73 

258 

3 

0.0588 

0.0136 

73 

259 

0 

0.0198 

0. 0154 

73 

259 

1 

0.0079 

0.0125 

73 

259 

2 

0.0060 

0.0058 

73 

259 

3 

0.0059 

0.0033 

73 

260 

0 

C. 0037 

0. 0046 

73 

260 

1 

0.0033 

0. 0034 

73 

260 

2 

0.. 0030 

0.0064 

73 

260 

3 

0.0050 

0. 0066 

73 

267 

1 

0.0467 

0. 0086 

73 

267 

2 

0.0289 

0.C126 

73 

267 

3 

C. 0426 

0.0436 

73 

268 

0 

0.0484 

0.0350 

73 

268 

1 

0.0659 

0.0366 

73 

268 

2 

0.0413 

0.0380 

73 

268 

3 

0.0065 

0.0417 

73 

273 

0 

0.0587 

0.0377 

73 

273 

l 

0.0226 

0.0376 

73 

273 

2 

0.0223 

0.0407 

73 

280 

0 

0.1117 

0.0399 

73 

200 

1 

0.0767 

0.0342 

73 

260 

2 

C.C499 

0.0292 

73 

280 

3 

0.0362 

0. 0293 

73 

281 

0 

0.0219 

0.0271 

73 

281 

1 

0.0177 

0.0255 

73 

281 

2 

0.0157 

0.0301 

73 

282 

0 

C. 0084 

0.0177 

73 

282 

1 

0.0070 

0. 0081 

73 

282 

2 

0.0091 

0.0049 

73 

232 

3 

0.0092 

0.0048 

73 

284 

3 

0.0052 

0. 0062 

73 

295 

0 

0.0849 


REPRODUCIBILITY OP ^ TT 

{ Mrf\\AJ PA op jp % 


TABLE D-l. Periods Used (continued) 


TIME 


Pin 

73 

295 

1 

0.0487 

73 

295 

2 

0. 0395 

73 

295 

3 

0.0411 

73 

296 

0 

0.0402 

73 

296 

1 

0.0304 

73 

296 

2 

0.0173 

73 

297 

0 

0. 0163 

73 

297 

l 

0.0038 

73 

297 

2 

0.0100 

73 

297 

3 

0.0336 

73 

298 

0 

0.0268 

73 

298 

3 

0.1091 

73 

299 

0 

0.0960 

73 

304 

3 

0.0136 

73 

305 

0 

0.0131 

73 

305 

1 

0.0109 

73 

305 

2 

0.0112 

73 

305 

3 

0.0080 

73 

317 

0 

0.0067 

73 

317 

1 

0.0123 

73 

317 

2 

0.0051 

73 

322 

0 

0.0032 

73 

322 

1 

0.0091 

73 

322 

2 

0.0073 

73 

322 

3 

0.0045 

73 

323 

0 

0.0032 


TIME 


PLO 

73 

323 

1 

0.0041 

73 

323 

2 

0.0073 

73 

323 

3 

0.0152 

73 

324 

0 

0.0137 

73 

329 

2 

0.2179 

73 

331 

1 

0.0525 

73 

331 

3 

0.0500 

73 

332 

0 

0.0450 

73 

332 

1 

0.0534 

73 

332 

2 

0.0611 

73 

332 

3 

0. 0648 

73 

333 

0 

0.0736 

73 

333 

l 

0. 0642 

73 

333 

2 

0.0711 

73 

335 

1 

0.0080 

73 

335 

2 

0. CO 74 

73 

335 

3 

0.0089 

73 

336 

1 

0.0134 

73 

336 

2 

0.C131 

73 

336 

3 

0.0154 

73 

342 

2 

0.1328 

73 

342 

3 

0.0804 

73 

343 

0 

0.0402 

73 

343 

1 

0.0336 

73 

343 

2 

0.0359 

73 

343 

3 

0.0264 


TIME P LO 


73 

344 

0 

0.0350 

73 

344 

1 

0.0683 

73 

344 

2 

0.0735 

73 

344 

3 

0.0968 

73 

345 

0 

0.1391 

7 3 

346 

1 

Q. 0072 

73 

346 

2 

0.0100 

73 

346 

3 

0.0077 

73 

347 

0 

0.0082 

73 

347 

1 

0.0125 

73 

347 

2 

C. 0147 

73 

347 

3 

0.0156 

73 

348 

0 

0.0120 

73 

348 

1 

0.0138 

73 

348 

2 

C.0140 

73 

356 

3 

0.0041 

73 

357 

0 

C. 0055 

73 

357 

3 

0.0108 

73 

358 

0 

0.0066 

73 

358 

1 

0.0071 

73 

350 

2 

0.0048 

73 

358 

3 

0.0031 

73 

359 

0 

0.0042 

73 

359 

1 

0.0049 

73 

360 

0 

0.0034 
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